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abstract 

space simulated , hetman, induced de.ctma.ions and Cesses in continuous fiber teiniotced com- 
posites were investigated with a mictomeehanics analysis. The investigation focused on 
ppmaiy areas. First, available explicit sessions ,ot Pricing me elective «*"**"* o, the,- 
^ expansion (CTE's) .0. a c<c^«e wae competed *h each enhet, and with a .inne Fernet. 
(FE) analysis, developed specifically fot this study. AnO^ cempatisons were made ,ot a w, e 
, ^off, be, fmarlx systems, and fxediCed vaeswae competed wim experimeaa daa. Al, a 
the analyses predicted neatly indentical values of the axial CTE, <t, , fot a given malen sy 

a, a me predict wete in good a^eenea w«h me expert data. Rosuits horn me 

c *vxr, nf o npneralized plane strain boundary value problem, 
FE analysis, and those from the solution of a generalized p 

othpr and with the experimental data for the transverse 
were in excellent agreement with each other, and witn 

CTE. », Less rigorous famutaions »a. in poor agreemea wim me expend data. 

The second area o, Ration footed on fhe detennir*ion of thetmaiiy indeed sttess fiefds in 
„ individual constituents. Stresses ^edicted ftcvn .he FE anaysis wete compared ,0 those P-e- 
diced from a closed-, rom solution ,0 the composite cylinder (CC) model, for two carbon fiber/epoxy 
composites. A gioba-loca formuiaiat. — iaminaed pfae theo. and FE anaysis. was 
used to determine the cesses in muMreCiona taminaes. Thermal, indeed damage ,mt,a,on 
^ ,s were also made. The type o, analysis (i.e. CC a FE) was shown to significantly attec 
me distributions and magnitudes of the predicted stresses. Thermally indeed marc stresses 
ie, eased In abaftae * w«h increase fiber vCume fracion. bu, wae nc a strong funebon 
a fiber properties. MultidireCiona [Os/ ± St lammaes had lager prediCed them y 
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matrix stresses than unidirectional ([0]) laminates, and these stresses increased with increasing 
lamination angle 6 . Thermally induced matrix failure predictions, using a maximum stress failure 
criterion based on the normal interfacial stress component and the measured transverse lamina 
strength, were in excellent agreement with experimental data. 
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1.0 Introduction 


Continuous graphite fiber reinforced composites are candidate materials for many space structures 
because of their capability for high stiffness, low coefficient of thermal expansion, and light weight. 
This combination of properties allows for the design of large stiff structures with minimal thermal 
distortions. Examples of applications currently under consideration include the large truss structure 
of NASA’s Space Station Freedom (SSF) W , Figure 1 , and the support structure and reflector panels 
for NASA’s Large Deployable Reflector (LDR) spacecraft^. The NASA Hubble Space Telescope, 
already built and scheduled for launch in 1990, makes extensive use of composites in the opt.cs 
metering and support structured (Figure 2). Although composite materials offer the potential for 
enhanced performance and significant weight savings, there are concerns about their long-term 
stability in the space environment. 

Materials in the space environment are exposed to ultraviolet and particulate radiation, atomic oxy- 
gen, vacuum, micrometeoroids, and large cyclic changes in temperature^ . The level of exposure 
to most of these parameters depends upon the orbit and protection systems on the spacecraft. For 
example, spacecraft ;.i low-earth-orbit (LEO), approximately 250 to 500 miles, are exposed to sig- 
rv- . , amounts of energetic atomic oxygen. In geostationary-earth-orbit (GEO), at approximately 
22, 000 miles, the levels of atomic oxygen are negligible, but there are significant amounts of par 
ticulate (electron and proton) radiation. Spacecraft are exposed to vacuum, ultraviolet radiation, 
and the threat of micrometeoroid impact in both orbits. 
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Figure 1 . NASA Space Station Freedom. 



Figure 2. NASA Hubble Space Telescope. 
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Temperature cycling also occurs In both orbits, however the frequency of the cycling does vary with 
me orbit, in LEO, spacecraft are exposed to approximately 16 cycles per day. A LEO spacecraft 
designed for a 30 year lifetime, such as SSF, will experience appmximatel, 1 75,000 thermal cycles, 
in GEO, the thermal cycling frequency is 1 per day. The maximum and minimum temperatures in 
a given cycle can be controlled by various thermal control/protection schemes, including active 
heating/cooling, and/or passive coatings and enclosures. The coating system pmposed for the 
buss structure tubes of SSF have a predicted "worst case" temperature cycle of ±150-F <’>. 

order to successfully design spacecraft structures with composite materials the effects of the 
space environment must be well understood. Most of the environmental parameters discussed 
above can affect the response and performance of composite materials during the design lifetime 
of the structure* 41 . Previous research has shown that large temperature changes and/or repeated 
thermal cycling can cause significant damage in composite materials and adversely affect their 
performance! 5 -' (Figure 3). The majority of past research has focused on analysis and testing of 
laminate response. However, the effects of constituent properties and microstructural characteris- 
tics (i.e. micromechanics) have no, been thoroughly investigamd for composites exposed ,0 the 

thermal environment of space. 

The objective of this research is to analyze the defomtations and stresses, induced by thermal 
loading, in continuous fiber reinforced composites using the properties and behavior of the con- 

stitutents (i.e. fiber, matrix, and intedace). The research will focus on two primary areas. First, 

available explicit expressions tor predicting effective coefficients of thermal expansion for the com- 
poslte will be con-pced with each other, and with a finite element analysis that was developed 
as pa,-, oi this research. Comparisons will be made for a wide range of fiber/matrix systems to 
determine the influence of the constituent pcopedies. All of the predictions will be compared with 
experimental data to assess the validity and shortcomings of the individual analyses. The remain- 
der of the research will focus on the determination of thermally Induced stress fields in the individual 
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[0 2 /90 2 ] s P75/934 Graphite/Epoxy 


Thermal 

strain 



Figure 3. Effects of thermal cycling on cross-ply graphlte/epoxy laminates. 


constituents. Stresses predicted from a finite element analysis will be compared to those predicted 
from a closed-form analysis available in the literature. The constituent properties and microstruc- 
ture that control the behavior will be identified and material modifications to improve the behavior 
will be suggested. A simple global/local formulation will be used to determine the influence of mul- 
tiple ply laminate constraints on the constituent stresses. Thermal stress calculations will be used 
to predict probable damage initiation locations, and the results will be compared to experimentally 
observed damage in polymer matrix composites. 


Four papers have been published by the author during the course of this research* 10-13 *. This 
report will contain details not included in these papers, as well as new unpublished results. 

Identification of commercial products in this report is provided to adequately describe the products 
and does not constitute official endorsement, expresses or implied, of such products or manufac- 
turers by NASA 
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2.0 Literature Review 


The micromechanics .iterators can he divided into two categories. The first de,s with the dev.- 
opment of explicit expressions for directly predicting the effective properties ot the composite from 
me propedies of the constituents. The second category covets analyses that provide detailed ,n 
formation about the displacements, strain, and stresses in me composite, rom which the e«ecr,e 
properties can he calculated. Extensive msearch fo, geneta. mtfliphase compose me.en,s has 
been conducted and published in both of these categories. This rev.ew w.ll be hmited to those pa 
pers dealing specifically with continuous fiber reinforced composite. Furthermore, this rev.ew w, 
address only those papers that deal specifically with predicting the effective coefficients of thermal 


expansion (CTE's), or present stress 


analyses that are applicable to thermally induced loading. 


2.1 Effective Coefficients of Thermal Expansion 


Numerous expressions have been derived. w„h voting degrees o, complexity, for predicting me 
m-piaie CTE's o, unidi.ectional comings fiber reWorced composite, A summary o, research 
„„ a,ng both the effective mechanical and thermal properties of composites was given by 
Cbamis and Sendee*, i<«>, where various Elyses were compamd with each other and experi- 
mental dam. Chamis and Sendeckyi divided the analysis methods into the following ma,or cate- 
gortes- strength o, mortals, se«-consis,ent models, vanatlonal or energy approaches, and exact 
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Figure 4. Unidirectional composite geometry and principal material coordinates. 

(i.e. using the concept of a periodic array of fibers and a repeating unit cell). Several theories have 
been proposed since this publication, and much more experimental data exists on current material 
systems. A bnef descnption of some of the more widely used current analyses will be given here. 

The geometry of the problem under consideration is shown in Figure 4. There are several basic 
assumptions that are common to all of the analyses to be presented and can be stated as follows: 

1. The fibers are circular in cross-section and infinitely long in the 1 -direction. 

2. All of the constituents exhibit linear elastic material behavior. 

3. The displacements are continuous across the fiber/matrix interface. 

4. The temperature distribution is uniform throughout the body, and the constituent material 
properties do not vary with temperature. 

ShaperyOS) derived expressions for the effective longitudinal and transverse CTE's (*, and <* 2 , 
respectively) of a unidirectional composite with isotropic fibers embedded in an isotropic matrix, 
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based on a variational energy appraoch. These expressions may be written as 


01 


EXV* + E m o m V m 


yi + E m V m 


(2.1.1) 


and 


a 2 = (1 + t' m )a m V m + (1 + - a 1t / 12 


( 2 . 12 ) 


The terms E, u, and V refer to Young’s modulus, Poisson’s ratio, and phase volume fraction, re- 
spectively. Terms without superscripts in equations (2.1.1) and (2.1.2) refer to effective composite 
properties. All other terms with superscripts of f and m, refer to fiber and matrix properties, re- 
spectively. This superscript convention is applied throughout the paper. The subscripts used in 
equations (2.1.1) and (2.2.2), and throughout this section, refer to the principal material coordinates 

of Figure 4. 

Unfortunately, most fibers such as graphite are not isotropic, but are orthotropic. Graphite fibers 
are often considered to be transversely isotropic (i.e. the plane of isotropy is perpendicular to the 
longitudinal axis of the fiber). A further discussion of the implications of fiber orthotropy will be 
given in Section 5.1. Strife and Prewo (16) have used a simple modification of equation (2.1.2) by 
replacing the isotropic fiber CTE, a \ , with the transverse fiber CTE, o f 2 , to account for transversely 
isotropic fibers. This modified expression is written as 


a 2 = (1 + r m )a m V m + (1 + *'{ 2 )a 2 V f - (2. 1 .3) 

However, the authors note that the modification is not based on any mathematical derivation. Com- 
parisons between experimental data for unidirectional Kevlar/epoxy and predicted values from 
equations (2.1.1) and (2.1.3) have been reported by Strife and Prewo< 16 >. A simple rule-of-mixtures 
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formula written as 


i/ 12 = i/{ 2 V» + ^ m V m (2.1.4) 

was used for the effective longitudinal Poisson’s ratio of the composite needed in equation (2.1.3). 
The agreement between experimental and predicted values of was very poor. The authors 
attributed the poor agreement to inaccurate input values of a’ . The predicted values of a 2 were 
within 15 percent of the measured values. 

Chamberlain* 1 ^ used the plane stress solution of a thick walled cylinder to derive simple expres- 
sions for a-i and a 2 of a transversely isotropic fiber embedded in a cylindrical isotropic matrix 
region. The expression for is identical to equation (2.1.1). The expression for a 2 takes the form 


m + M -°" l > Vl 

2 a „«"(F- -I +v m ) + (F + V*) + (E m /E < 1 )(1 -I^ 2 )(F-1 +V m ) 


(2.1.5) 


where F is a packing factor which accounts for fiber packing geometry, and is equal to 0.9069 and 
0.7854 for hexagonal and square packing geometries, respectively. 

Rogers et al* 18 > used equations (2.1.1) and (2.1.5) to predict the CTE’s of unidirectional graphite/ 
epoxy composites, and compared the predicted values with measured values. The predicted val- 
ues were within 2 percent of the measured vaules of a 2 , using the hexagonal packing factor. The 
predictions, in absolute terms, were within 0.03 x 10 _6 /°F of the measured values. Strife and 
Prewo* 16) also used equation (2.1.5) to predict a 2 for Kevlar/epoxy, and reported differences of 
about 13 percent from the measured values. 
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used a simple force-balance, o, strengtfvo,- — ^ 

for „ CTE's o. unidirectional composes cons^ing C transvers* iscsropic tide, in an -trope 

matrix. The expression for «i is again i 
written as 


, identical to equation (2.1.1). Hie expression lor o 2 can be 


q 2 = a^V* + a m V m (1 + V*n m E , 1 /Ei) 

where Ei is the longitudinal Young's modulus of the composite, given by the simple rule-of-mixtures 
formula written as 


E 1 = E’y + E m v r 


(2.1.7) 


:’S 


No comparisons between measured and predicted values using equation (2.1 .6) could be lound 
in the literature. However, a comparison between equation (2.1.6) and a Unite element anay 

showed very poor agreement^ . 

Rosen and Hashin (21 1 extended the work ot Levin' 22 ! to derive expressions (or the effective CTE 
o, multiphase compost. Hashing summed this approver unidirediona, liber reinforced 
composites w«h transveray isotropic constituent This anrSysis is the most genera, ol the ana,- 
yse s discussed thus tar. and relates the volume average stresses and strains in a characteristic 
volume element to the effective properties of that element. The derivation leads to comp 
loose, expressions for and and requires three of the effective elastic mechanical properties 
(longitudinal modulus and Poisson, -aio, ami transverse « moduli*, o, the composite in the 
calculations. In contra* the Shapery an*, as requires on* the ertectNe longrtudina, Poissons 
ratio, die Chamis analysis requires only me effect longitudinal modulus, and me Chambertam 
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analysis requires none of the effective elastic mechanical properties. The general form of the results 
may be written in tensor notation as< 23 > 


<*ij - «ij + (of - «ij m, )Pklrs(S rs jj - S rsi j) (2.1.8) 

P klrs(Sreij ~ S rij) = •klij (2.1.9) 

where S rsj j are the elastic compliances and l kM j is a fourth rank symmetric unit tensor. Superscripts 
(0 and (m) refer to fiber and matrix, respectively, are enclosed in parentheses to distinguish them 
from tensor indices. The terms with an overbar and hat refer to effective and volume average 
composite properties, respectively. Composite volume average properties are obtained by 


p = v f p f + V m p m 


( 2 . 1 . 10 ) 


where p is the property of interest. 


For transversely isotropic constituents, equation (2.1.8) simplifies to 


Q 1 01 + ( a k! ~ Q kl *) P klrs(S rs i1 - S rs11 ) (2.1.11a) 

Q 2 = a 2 + (»S - a&' 4 )P|dm(Sia22 ~ S rs2 2) (2.1.11b) 

The overbar has been dropped from qi and a 2 to be consistent with equations (2.1.1)-(2.1.7). 
Equations (21.11a) and (2.1.11b) are solved by determining P Wr# from equation (2.1.9). The 
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components of S are obtained from expressions for the effective mechanical properties given by 
Hashin* 23 ' Equations (2. 1 . 1 1 a) and (2. 1 . 1 1 b) were expanded and simplified to a form more suit- 
able for routine computations in Appendix A. To the author’s best knowledge, this is the first time 
that a fully expanded version of equations (2.1.11a) and (2.1.11b) has appeared in the open liter- 
ature. No comparisons between measured and predicted values using equations (2.1.11a) and 
(2. 1 .1 1 b) could be found in the literature. This may be partly due to the complexity of the expres- 
sions relative to some of the other previously discussed derivations. 

Exact methods employing the assumption of a regular periodic array of fibers, from which a re- 
peating unit cell may be extracted, have also been used to predict CTE’s. A discussion of these 
methods will be given in the section on stress analysis formulations. 

As described above, numerous expressions have been derived, with varying degrees of complex- 
ity, for explicitly computing the effective CTE’s of unidirectional composites, based on the elastic 
properties and CTE’s of the constituents. However, a systematic comparison of these different 
formulations with each other, and with a broad range of experimental data has not been reported. 
This type of comparative study is essential for determining the applicability of each formulation, 
and will be one of the major thrusts of this dissertation. 

2.2 Thermal Stress Analysis 

Tne prediction of thermally induced damage and comparisons with experimental data have, for the 
most part, been limits' 1 ;^the macro or laminate level^ 9-10 ). However, numerous papers have been 
published on predicting the mechanical stress-strain behavior of unidirectional composites using 
micromechanics analyses. The three most common approaches fit into the three categories listed 
in Section 2.1 as exact, self consistent, and strength of materials formulations. Representative 
papers from all three categories are discussed below. 
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Exact formulations usually assume that the fibers occur in a regular periodic array. This allows the 
analyses to be performed on a repeating unit cell rather than the whole array of fibers, thus greatly 
reducing the size and complexity of the problem. The influence of adjacent fibers is incorporated 
through the boundary conditions on the unit cell. Many of the exact formulations have employed 
the finite element method. Adams et. al.* 24 ' 27 * have published extensively in this area using a 
generalized plane strain formulation, with constant strain triangular (CS1) elements, to model cir- 
cular fibers in a square array surrounded by matrix material. This analysis was used to examine the 
effects of matrix material nonlinearity and temperature dependence on graphite/epoxy composite 
stress- strain behavior for various load conditions. All of the thermal load conditions were for room 
temperature and above. Foye< 28 > and Dvorak et. al.* 29 * have also used finite element analyses with 
CST elements to predict nonlinear stress-strain behavior of unidirectional composites subjected to 
thermomechanical loading. These papers were primarily aimed at predicting matrix yielding in 
metal matrix composites subjected to mechanical loading. 

Aboudi* 30 - 32 * has also used an exact formulation. The Aboudi model assumes a triply periodic 
array of parallelepiped fibers embedded in an infinite matrix region. The representative cell for this 
model is a single parallelepiped fiber subcell surounded by three parallelepiped matrix subcells. 
This general geometry can be used to model particulate reinforcements, short fibers, continuous 
fibers, and 3D woven fabrics. A first order displacement expansion is employed in each subcell, 
which leads to average constant stresses and strains in each subcell. Imposition of continuity 
of tractions and displacements across the boundaries of the individual subcells leads to closed 
form expressions relating the applied stress to the average stresses and strains in each subcell. 
This model has had very good success in predicting the global effective properties and stress- 
strain behavior of metal matrix composites* 33 * . However, the assumption of parallelepiped fibers, 
and constant subcell stresses and strains, raises concerns about the accuracy of this model for 
predicting the local stress fields in the vicinity of the fibers. 
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Several stress analyses have been formulated using variations of the self consist 
of the self-consistent models assume that the real composite can be replaced by a representative 
cell embedded in a homogeneous medium whose properties are equivalent to the effective proper- 
lies of the composite. TN*e modNs neglect certain microstrucaaa. de«s, sr*h as the influence 
o, adiacen, *bers. «. are usual, best suited for prediction N ^oba, behavior. One of me most 
wide, used self -consistent models is me compose cylirxfer assemblage (CCA, -emulation o, 
Hashing The represent cel, in this formulation is a compose cylinder (i.e circular fiber em- 
bedded in a cylindrical region o, matrix) and can be so,ed as an etestierty boundary value problem^ 

A thorough description of this formulation will be given in Section 3.3.2. Avery and Herakovrch 
have used this method to examine the effects of fiber anisotropy on the thermomechanical stresses 

in unidirectional composites. Mlkraa «. Taya<*> have also used , Ns method to evamlne the ef- 
fete of fiber coalings on thermomechanicN saesses. Bom of mese ana,ses wem linear elastic 
with temperature IndeperrteN maeri* progenies, amt d« no, consider thermN loads below room 
temperature or the initiation of thermally induced damage. 

Another variation o, the sek-corsisten, mod* is me 'sNid mfcaae- —ion of Min<">. This is 
a pu* stress elasticplastic ana,sis W»h un«orm saesses in me fber and maaix phases. This 
andysis was used to study the eiaaicplastic mechanic* msponse o, metal matnx composites. 
Dvorak et al* 38 ' have derived a ’’vanishing fiber diameter” formulation which is also a variat 
the self-consistent model. The composite is modeled as acontinuum reinforced by cylindrical fibers 
of vanishing, small diametem tha, occupy a finite volume fraction of the composSe. This leads to 
uniform local stress c..d strain fields in me consents, bu, great, simpLies the oomputa.- o, 
ms vi composite stress-strain behavior. Tl» atshom do aate that this formuWon neglecN 
certain INeracSons behveen f*ases in the averse direcSon, and therefore has 
accuracy for predicting composite response in that direction This ana,sis was also used to study 

me elastic-plastic mechanical behavior of metal matrix composites. 
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Chamis (39) has used a strength of materials approach to derive simple expressions for fiber and 
matnx stresses using constituent properties, effective unidirectional lamina properties, and applied 
lamina stresses. Fiber stresses are assumed to be uniform, and there are two regions of uniform 
matrix stresses. This formulation is consistent with the formulation given by Charm's* 1 9 > for predict- 
ing effective properties. 

Analyses that fall under the category of exact, as described at the beginning of this section, give 
more accurate descriptions of the local variations in the stress and strain fields of the constituents 
by virtue of their formulation. Analyses based on variations of the self-consistent model have been 
shown to be well suited for predicting global properties and response, but lack the microstructural 
details to predict localized phenomenon. If the unidirectional composite response is needed as 
input to a laminate model, then a global response might be adequate and even desirable due 
to its usually simpler computational scheme. However, accurate predictions of thermally induced 
damage require detailed information about local variations in the stress field, and thus require 
analyses that fall within the exact category. None of the micromechanics studies found in the 
literature, and described above, investigated the initiation of thermally induced damage in polymer 
matrix composites. An analytical investigation of this problem will be one of the main thrusts of this 
study. The results of an exact formulation based on a finite element analysis will be compared to a 
self-consistent scheme based on the composite cylinder model. 
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3.0 Theoretical Development 


3.1 Approach 

The theoretical thrust of this dissertation is divided into two areas. The first deals with the prediction 
of effective CTE’s of unidirectional composites. Four formulations for explicitly predicting effective 
CTE’s were presented in Section 2.1, and a discussion of these formulations will not be repeated 
in this chapter. Those expressions assumed temperature independent constituent properties. A 
method for correctly incorporating the effects of temperature dependent constituent properties will 
be developed and presented in this chapter. 

The second area of theoretical investigation deals with the prediction of thermally induced stresses 
in the constituents of unidirectional composites. Two formulations will be presented in this chapter; 
a finite element analysis, and an elasticity solution of the composite cylinder problem discussed 
in Section 2.2. The*" two formulations will be presented for temperature independent constituent 
rr js first. A general method for incorporating the effects of temperature dependent con- 
stituents will then be presented. A brief discussion on the use of classical laminated plate theory 
(LPT) for predicting laminate thermal response, and a simple global/local formulation combining 
LPT with micromechanics will also be presented. 
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The general problem geometry under consideration was shown in Figure 4. The assumptions that 
are common to all of the theoretical work of this study are listed below. 

1. The fibers are circular in cross-section and infinitely long in the 1 -direction. 

2. All of the constituents are linear elastic. 

3. The displacements are continuous across the fiber/matrix interface. 

4. The temperature distribution is uniform throughout the body. 

These assumptions are the same as those used in the derivations for predicting effective CTE’s 
discussed in Section 2.1 , with the exception that the constituent properties may now be functions 
of temperature. 

3.2 Effective Coefficients of Thermal Expansion 


The expressions for predicting effective CTE’s presented in Section 2.1 were derived assuming 
temperature independent constituent properties. If the constituent properties vary with temperature 
the expressions are no longer valid, and must be modified. The appropriate modification consists 
of replacing all of the r*’s (effective, fiber, and matrix) with thermal strains, e, given by 

/r, 

f(T) = / a(r)dr (3.2.1) 

-'Tsft 


where T sft and T, are the stress free temperature and analysis temperature, respectively. The 
Tsft is usually assumed to be the temperature at which the composite is fabricated, and as the 
name implies, is the temperature where the composite and its constituents are all stress free. 
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Conversely, by definition a is given by 


q(T ) _ dT 


(3.2.2) 


t ic nhtained from the derivative 

Therefore, the instantaneous effective GTE of the compose at T, - obtained 


of the expression for the 


effective thermal strain of the composite evaluated at T, . 


The above procedure is best illustrated by 


considering a specific example. The expression for the 


effective 


ive longitudinal , harms, strain is -mm equaion 2.1.1 and may be wrtten as 


+ E m e m V m 
f1 = £» v< + E m V m 


where all ot the constituent propenies may now be lunctions ot temperature. 

The effective iongffudina CTE is now computed horn ecffiaffon (3-2.2, and is wdtten as 


(3.2.3) 


«1 = 


d ( Eyy +£ m t rn v m \ 

dT l + E m V m "" ) 


(3.2.4) 


Now, if the moduli, E*, and E m , are 
the form of equation (2.1.1), written as 


not lunctions of temperaure, men equation (3.2.4) simplifies to 


+ E m a m V m 

a 1 = yl + Em\/m 


where a\ and are the CTE's of the constituents 


(3.2.5) 

evaluated at the temperature of interest, Ti . 


However, if th~ c*,. .stituent moduli are functions 


of temperature, then equation (3.2.4) becomes 

dE m \ 


a 1 


jpf dE m \ 

_ __j ( E y a \ + E m v m « m + ) 

- £f vf + E m V m V ' 


2 r 


E f\/f + E m V m 


x l f dEt . <m dE m ^ 

(eWM+E-V^^-^+V’ -fir 


(32.6) 
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Equaion (3.2.6) IS raher lengthy, and its evaluaion a a particular temperature may be vety cum- 
betaome. This Is especiaiy hue if the functiona form a the constituent property vartalon with 
temperaure is not Known The development o, explicit expressions for the temperature depen- 
dent transverse CTE becomes even more lengthy and complicated. 

An aternae apptoach to determining explicit expressions for the temperaunr dependent effective 
CTE-S. as outlined above, is to lira compute effective therma strains, and then numericaiy evalu- 

aemederivadveae q uation(3.2.2). T heeffeaive,hermaartKnsaecompu,edbysubai,utin 9 me 

conaituent thermal arains (defined by equaion (3.2. 1» fo, the conakuen, CTE's into the expres- 

aorts presented in Section 2.1. An example of this subaitution was shown in equaion (3.2,3). 

This computaion requires tha the conaituent elaaic properties be known only a the analysis 

temperaure. The functional form of their vaiaion with temperature is not required. However, the 

pah dependence of the constituent CTE’s with temperaure is required for evaluaion of equaion 
(3.2.1). 

The numerical approach used in this study for evaluaing equation (3.2.2) was to fit a second order 
interpolaing polynomia to three discrete effective therma arain values. The three temperatures 
selected wae the analysis temperature, T, , and temperatures T, ± 5-F . The expression for the 
derivative of the effective thermal strain evauated a T, may then be written as«°) 

= df = To ^ + 5 ) “ - 5 )) (3.2.7) 

The above expression was usedfor both the longitudinal and transverse CTE'saeach temperature 
of interest. 
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3.3 Thermal Stress Analysis 


3.3.1 Finite Element Model 

Afinite element (FE) model was used to determine the thermally induced stresses and deformations 
in unidirectional composites. The FE analysis of the unidirectional composite geometry, shown 
earlier in Figure 4, was greatly simplified by assuming the fibers are arranged in some type of regular 
and periodic array. This allows the analysis to be performed on a representative unit cell rather 
than the entire composite, and greatly reduces the size and complexity of the analysis. Additional 
simplifications to the unit cell can be obtained by utilizing the symmetry of its geometry and applied 
loads. The two fiber array geometries assumed for this study were square and hexagonal. The 
simplified unit cells, boundary conditions for a uniform thermal load, and coordinate systems for 
these two array geometries are shown in Figure 5. The (1,2,3) principal material coordinates of 
Figure 4 are coincident with the (x,y,z) cartesian coordinates of Figure 5. 

A condition of generalized plane strain was assumed to exist in the composite. The displacement 
field for this condition may be written as 

u = U(y,z) + x 

v = V(y,z) (3.3.1. 1) 

w = W(y,z) 

where U, V, and W a;e unknown functions of the spatial coordinates y and z, and is a uniform 
strain (i.e. constant in the yz plane) in the x coordinate direction. The term may be a known 
applied strain or treated as an unknown for a uniformly applied force, FJ in the x direction. For 
thermal loading only, the uniform axial applied force is zero and is the unkown thermal strain in 

the x direction. 
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P = 30 deg for hexagonal array, 45 deg for square array 

Figure 5. Finite element micromechanics model geometry and boundary conditions. 


A linear elastic displacement formulation was used to solve for the above unknowns, using eight 
node isoparametric quadratic elements. The FE mesh used in this study is shown in Figure 6. 
This final mesh geometry was derived from a mesh refinement study to determine the minimum 
number of degrees of freedom (DOF’s) that would provide convergence of the solution for the 
displacements. Details of displacement based FE formulations are found readily in the literature^ 41 ). 
A description of the formulation is given in Appendix B. This description covers the solution for the 
nodal displacements as well as the computation of the element stresses and strains from these 
displacements. 

Several unique features, specifically required for this study, were incorporated into the formulation 
presented in Appendix B. First, the inclusion of an unkown is carried throughout the formulation, 

and its effect on the global system of equations is highlighted. A detailed treatment of this type 
could not be found in the literature. A method to account for material cylindrical orthotropy is also 
included in the formulation. This type of orthotropy caused the material stiffness matrix to vary 
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8-node quadratic isoparametric elements 
158 nodes, 45 elements 


Fiber/matrix interface 



Figure 6. Mlcromechanlcs finite element mesh. 

^ a given elemem. The Nation was inciuded in me formuWon W ev*u*,ng me stress 
matrix a. each of the integration points used in me numerical integration scheme tor an element. 

F ^, aspect solution algorimrn was damped, or e« simringthe ^oha, system otiine, 

equations for the case of an unknown r J. This was necessary because me inclusion of an unkown 
destroys me handed of ft. glob* system * er^tens. The specif sloped a,gor # hm 
allowed the h^ded nature o, the system to he retained. ftu. greedy basing sofuion ehicmncy. 

The formuiation descdhed ahove and in Appends B was imp.emen.ed in a compeer program 
called FECAP (Finite Element Composite Analysis Program) .and was written to run interactively on 
a desktop microcomputer* 111 . FECAP was written in Hewleh Packard BASIC 3.0, and can run on 
any wrett Packard 9000 Series microcomputer with a BASIC 3.0 or compatible operating system. 

Nodal coordinate and element connectivity date are read from a user generated file. A very simple 

mesh generation program was written to produce me data files used in mis study. The output horn 

FECAP (i.e displacements, stresses, and drains, was wrWen to a fife for Esther po*.proces*ng 
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and plotting. 


3.3.2 Composite Cylinder Model 


The other analysis used for predicting thermally induced stress fields was the solution to the com- 
posite cylinder (CC) boundary-value problem. The CC model, which consists of a circular fiber 
embedded in a cylindrical region of matrix, is shown in Figure 7, along with the coordinate sys- 
tern and boundaty conditions. The use of the CC model for providing the solution to a composite 
nsisting of many fibers (Figure 4) is based on the composite cylinder assemblage (CCA) for- 
mulation of Hashing. This formulation, a variation of the sell-consistent approach, assumes that 
a composite can be modeled as a collection of composite cylinders of vaiying size, all of which 
have the same ratio of fiber radius lo mairix radius (i.e fiber volume fraction). A single composite 
cyfinder becomes representative of the entire composite by esquiring that, for a given load state, 
the stored strain energy in this single composite cylinder is equal to the strain enenjy in a cylinder 
of homogeneous material with the same "effective- thermoelastic progenies. Consequently, the 

stresses and strains in a single composite cylinder represent only the average stresses and strains 
in the “real” composite. 


The CC model is most easily solved using the cylindrical coordinate system shown in Figure 7 . The 
terms u.,andw in the following formulation refer to displacements in the x, 9 . and r directions, 
respectively. Due to the axial symmetry of the geometty and load, for a uniform change in tem- 
perature, the V displacement will be zero and the solution will be axisymmetric (i.e independent of 
the 9 coordinate). Conditions of generalized plane strain may also be assumed for this problem, 
leading to a displacement field which has the functional form of 

u = u(r) + f ° - x 


w = w(r) 
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Figure 7. Composite cylinder model geometry and boundary conditions. 

where u(r) and w(r) represent the functional form of the u and v displacements, respectively, and 
t o i S an unknown uniform strain in the x direction. The strain-displacement relations can be written 


<9u 

d w 

W 


dx' 

ff " It' 

€ ' = 7 ’ 

(3.32.2) 


X 

II 

O 

du 


7f» = °> 

7xr " Tr 



The nonvanishing equilibrium equations for this problem are 

— + -(<r, - c,) = 0 P-3 23) 

dr r 

& t xt ^xr _ q (3.3. 2.4) 

~dr~ r 

Substituting the straiivdisplacement relations (equation (3.3.2 2)) into the appropriate material con- 
stitutive relations, and then substituting the resulting relations into equation (3.3.2.3) yields the 
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governing differential equation for the w displacement written as 


~ Cgf 2 - -(C #x - Cn<)e® + — (C r j - C # j)ojZiT (3.3.2 5) 

where the Cy and aj are the stiffness and CTE coefficients, respectively, and are defined in Ap- 
pendix D. The repeated subscript i in equation (3.3.2.S) is summed over x, r, and 0. The general 
solution to equation (3.3.2.S) has been given by Cohen and Hyer< 42 > for an orthotropic tube, and 
by Avery and Herakovich* 31 * for a fiber/matrix composite, and is repeated in Appendix D for com- 
pleteness. 


3 2 w 1 <9w 

3r 2 + r dr 


The solution of equation (3.3.2.4) may be written down by inspection as 


K 


r xr — 


(3.3.26) 


where K is a constant that can be determined form the boundary condition on r xr at the outer 
surface of the cylinder. For thermal loading only, the outer surface is stress free and therefore K 
must equal zero. This implies that there are no shear stresses or strains for the case of thermal 
loading only. Also, the functional form of u given in equation (3.3.2. 1) may be simplified to 


u = e x • x (3.3.2 7) 

where the displacement u is no longer a function of r. 


The solution to equation (3.3.2.5) and the subsequent expressions for the stresses in the fiber and 
matrix presented in Appendix D were coded into a computer program for ease of computation. 
This program was also written to run on a Hewlett Packard 9000 Series microcomputer. 
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3.3.3 Global/Local Model 


The boundaiy conditions tor the FE model shown in Figure 5 are tor a unidirectional composite 
with no externally applied loads except tor a Uniterm change in temperature. These boundary 
conditions must be modified to model an array of fibers and matrix in a laminated composite, in 
order to accountfor the constraints imposed on the individual laminae. These constraints arise from 
me mismatch in CTE (., and a 2 ) between plies with dilterent fiber orientations when the laminate 
is subiected to change in temperature. The therxy governing mis behavior in thin laminates, where 
conditions ot plane stress may be assumed, is known as laminated plate theoty (LPT) and is found 
readily in the literature <«>. In a symmetric laminate subjected to a uniform change in temperature 
the strains in each ply are equal, and may be written as 


1 1 1 - tAJ- 1 (£[6k W a] ^ < 3 ' 3 3 

where <„}„ is the CTE mask ot the klh ply, t k is the thickness of the k'th ply, and jAj and [Qj, 
are the stiffness matrices of the laminate and k'th ply, respectively. The summation is over the total 
number of plies, n. Definitions of the above .etms can be found in the LPT reference cited eadier. 

The strains in equation (3.3.3. 1) are referenced tome (x,y,z) laminate coordinate system shown in 
Figure 8, and must be transformed to the (1 ,2.3) principal material coordinate system (Figure 4) for 
use with the FE analysis. This is accomplished through a simple transformation equation written 


as 



cos 2 # 

sin 2 # 

-2cos#sin# 


sin 2 # 

cos 2 # 

2cos#sin# 


cos#sin# 
-cos#sin# 
cos 2 # — sin 2 # k 



where # is the angle measured from the x-axis to the 1 -axis for the k’th ply. 


(3.3.32) 
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The through-the-thickness (3 or z direction) thermal strain of a given ply is also constrained in a 
multidirectional laminate. The expression for this strain is derived directly from Hooke’s Law, and 
is written for the k’th ply as< 44) 


f3 = a 2 AT - <r, - <r 2 (3.3 3.3) 

The above equation assumes that the lamina or ply properties are transversely isotropic (2-3 plane 
of isotropy). The terms ^ 12 and i/ 23 are the lamina Poisson's ratios, Ei and E 2 are the lamina 
moduli, and a 2 is the unconstrained lamina CTE. The terms a 1 and <r 2 are the thermally induced 
lamina stresses in the principal material coordinates. These stresses are determined from LPT, and 
are written as 


| Z \ =[T] k [Q] k ({ f °}-{«}k^T) (33.3.4) 

l r i2 J k 

where {c 0 } is defined by equation (3.3.3.1), and [T] k is a transformation matrix given by 


(T] k = 


cos 2 # sin 2 # 2cos#sin# 
sin 2 # cos 2 # -2cos#sin# 
-cos#sin# cos#sin# cos 2 # - sin 2 #] k 


(3.3.35) 


The method for imposing the laminate induced constraints on a fiber/matrix unit cell is shown 
schematically in Figure 8. The thermally induced laminate strains, e j , f 2 , and 7 12 , and the through- 
the-thickness lamina strain, f 3 , were determined from equations (3.3.3.1) through (3.3.3.S), for a 
given AT, and ply. These strains were then used as displacement boundary conditions on the FE 
model of a single fiber and surrounding matrix exposed to the same AT. The FE model shown 
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Multidirectional laminate 
ply strains 

determined from LPT 
(e x ,e y ,e z ,Y X y) 


t z 

£ 

Matrix ^ 


Fiber^^^ 

— 


Idealized 
fiber 
array 


Boundary conditions 
from ply strains 


* Finite element model 

Figure 8. Finite element model for global/local formulation. 

in Figure 5 could no. be used for this formulation because of a dHference In symmelry conditions, 
instead, the -quarter-symme,*- mode, (i.e. modefing o' a quaner of a fiber) shown in Figure 8 
was used. The axial drain condition was imposed by sening me uniform applied axial strain, r? 
o, equafion (3.3.1.1) irom ihe FE formulation, aqua, to The transverse strain condition was 
imposed by setting the V displacement along the y=1 edge equal to < 2 . The V displacement along 
y=0 was set equal to zero. The shear strain condition was imposed by setting the U displacement 
^ y=0 to zero, and ihe U dispiacemed dong y-1 to y,, Findly. the W dispiacemen, dong 
, = n was set equal to zero, and the W displacement along z-1 was set equal to , 3 . 


3 3 4 T~ m perature Dependent Constituents 

The CC and FE models described above were presented for the case of temperature independent 
constituent propenies. However, both modds wem implemented w«h the capabilrty to dlow for 
temperdure dependent constituent mdend properties. The formulations were modified by follow- 
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ing the procedure described in Section 3.2. Terms containing a AT were replaced with thermal 
strains given by 

/ Tl 

e(T) = / a(r)dr (3.3.4. 1) 

•'T SFT 

which is identical to equation (3.2.1). The stresses in the constituents may be written in generic 
form as 

/•t, 

= H(Tl) / G(r)dr (3.3.4.2) 

■'Tsft 

where HO - ! ) is some function of the constituent elastic properties evaluated at T, , and G(r) is some 
function of the constituent CTE’s. 

This formulation is refered to in the literature as a “total strain theory"* 45 ). As discussed in Section 
3.2, the elastic moduli and Poisson’s ratios are only required at the analysis temperature. Their path 
dependence with temperature is not required. This formulation is only valid for elastic constituent 
material properties (i.e. properties not a function of stress). When the properties are inelastic, an 
incremental approach must be used. 
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4.0 Experimental Methods 


A brief description of the experimental method used to determine the thermal strain behavior of 
composites will be presented first. This will be followed by a discussion of the method used to 
observe and characterize thermally induced damage in a high modulus graphite/epoxy composite. 

4.1 Measurement of Thermal Strains 

Thermal strain data were determined for a number of different continuous fiber reinforced unidirec- 
tional composite systems. These data will be presented in Section 5.2. All of the measurements 
were made using a laser interferometric dilatometer at NASA-LaRC. This equipment has a strain 
resolution of approximately 1 xl (T 6 and a temperature range of - 250 0 F to +300 0 F . The specimen 
length is approximately 3 inches, and the strain is determined over the full length of the specimen. 
All of the resin matrix composites were vacuum dried to a constant mass before testing, and all 
of the specimens were tested in a dry N 2 environment. Details of the experimental technique are 
given in the literature* 46 ' . 

Vaiues of Ql and cv 2 were obtained by fitting a least-squares polynomial (usually 3rd order) to 
the thermal strain-temperature curves and then evaluating the derivative of this polynomial at the 
desired temperature. A typical axial thermal strain-temperature curve for unidirectional T300/5208, 
along with the least-squares polynomial fit, is shown in Figure 9. 
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Figure 9. Axial thermal strain curve for unidirectional T300/5208. 

4.2 Damage Observation and Characterization 

The formation and growth of thermally induced damage was observed and characterized in a high 
modulus P75/934 graphite/epoxy composite system. A description of the constituent properties 
of this particular system will be given in Section 5.1. Unidirectional and cross-ply panels were 
fabricated from prepreg tape according to the manufacturer's standard cure cycle. The maximum 
temperature in this cure cycle was 350 °F , which was the assumed stress-free-temperature in the 
thermal stress analysis to be presented in Section 5.3. The panels were ultrasonically C-scanned 
after fabrication, and exhibited r>o significant voids or delaminations. Both panels had an average 
fiber volume fraction of 0.50. 


Specimens from unidirectional and cross-ply laminates were observed after fabrication, and af- 
ter various numbers of thermal cycles between ±250 °F , for the presence of thermally induced 
damage. The observations consisted of examining the polished edges of the specimens with an 


Experimental Methods 


30 


optical microscope a. magnetons C <00* The spoken ^mansions were 3» by 1 M * 
0 0< m inches. M o. the specimen were deed pier to obsecabon and were kept ary «tro^h- 
out ft. thermal cycling process. The edects o. therma, cyciing intkKed damage on the module* 
and GTE ot these specimens, as well as the details o. the experimental procedure, have been re- 
ported elsewhere 1 * 0 ’ . The present study will focus on die initiation and morphology of the damage 


observed in these specimens. 
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5.0 Results and Discussion 


The discussion of the analytical results will be divided into four main sections. The first section will 
discuss the constituent property data that was used as input for all of the analytical predictions. 
The second section will discuss the various methods for predicting effective CTE’s of unidirectional 
composites, and their comparisons with experimental data The third section will describe the 
results of the thermal stress predictions, and the final section will discuss the formation of thermally 
induced damage. 

5.1 Constituent Properties 

As discussed in Section 2.1, carbon fibers are not isotropic, but are commonly assumed to 
be transversely isotropic with the plane of isotropy being the 2-3 plane of Figure 4. How- 
ever, studies* 47-51 ) of carbon fiber morphology suggest that the fibers may posses a combi- 
nation of a cylindrical^ orthotropic (i.e different properties in the x, r, and 0 directions of Fig- 
ure 7) sheath near the outer surface of the fiber, and a transversely isotropic core as shown in 
Figure 10. The type of cylindrical orthotropy (i.e. radial or circumferential) depends upon the fiber 
fabrication process. Carbon fibers made with a PAN (polyacronitrile) precursor typically have a 
circumferentially orthotropic sheath and those made from a pitch precursor typically have a radi- 
ally orthotropic sheath. A thorough investigation of the effects of fiber morphology on composite 
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Transversely Isotropic (Tl) 



Circumferentially Orthotropic (CO) 



Radially Orthotropic (RO) 


Figure 10. Carbon fiber morphology. 


eflecfce ,*ope„ies «l stresses has been given by Knob**. Specie resu»s hem the, study will 
be discussed in subsequent sections. 

The majority of results for the present study were generated with the assumption of transversely 
isotropic fibers. This reduces the number of independent elastic constants for the fiber to five, 

written as 

E \ . E 2 = E 3’ "12 = "13. "23. G 12 = G< 13 
with the transverse shear modulus related to the other properties by 

g£ 3 = 4/2(1 + 4 3 ) 


(5.1.1) 


( 5 . 1 . 2 ) 


All of the matrix 


materials were assumed to behave 


isotropically With two independent elastic 


constants written as 


( 5 . 1 . 3 ) 
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and the shear modulus defined by 


G m = E m /2(1 + ^ 


(5.1.4) 


E, i/. and G in the above expressions represent the Young’s moduli, Poisson’s ratios, and shear 
moduli, respectively. The subscripts refer to the material coordinates of Figure 4, and the super- 
scripts f and m refer to fiber and matrix, respectively. In addition to the elastic constants, there are 
two independent CTE's for a transversely isotropic fiber, a f t and a f 2 , and one for an isotropic matrix, 

Or 171 . 

Fiber property data, especially as a function of temperature, are very difficult to find in the literature 
due to the extreme difficulty of generating these values experimentally . The carbon fiber properties 
used in the present study were assumed to be temperature independent. This assumption is 
thought to be valid due to the relatively small temperature range used for the analyses (-150°F 
to +250 °F ) compared to the large useful operating temperature range of carbon fibers, and the 
large distance of the analysis temperature range from the fabrication temperatures of the fibers. 
The axial CTE of PAN based carbon fibers has been shown to be independent of temperature over 
this range* 53 ). The elastic constants and CTE’s for the various carbon fibers used in the analyses 
are given in Table 1. The T300, C6000, and HMS fibers are PAN based, and are manufactured 
by Amoco, BASF, and Hercules, respectively. P75 and PI 00 are high modulus pitch based fibers 
manufactured by Amoco. The data in Table 1 come from various literature sources, including 
both research papers and manufacturers product data sheets* 54-57 ). Experimentally measured 
values from the literature were used when available. However, the transverse fiber properties 
represent values from the literature that were back-calculated from composite properties using 
micromechanics theories for predicting effective composite properties* 23 ). 
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a 2 


El Ez G '\ 3«, 1/12 (io ~j/*g) (i°~V 0F > 

560 - 

I S « a SS ss 



(1) G23 = Ez/2(1 + ^23) 

(2) C6000 assumed to have same prope les 


5.60 

3.80 

3.80 

3.80 


A wide range of matrix materials was used (of the effective CTH predictions. The properties of 

were assumed to be temperature independent asswwed to be temperature 

taken Iron, the literature***. The polyimide propel 

— — — — — - 7^1 

F „»c» „ 

used for the majority of the stress analysis predictions. The other epoxies in Table 2 armc , 
d L CE339) «re * assumed to have the s*ne — ^ 

934 matrix properties were assumed to be temperssure dependent. Experimental values ol E - 

934 P H • Qn hu Fo v et al < 59) . The elastic properties 

/ot\ onH 4 - 950 °F were given by rox ex.<«. 

„rr. at room temperature (RT) and +250 i- , y (60) 

3 , were hack-caicuiaied to, «. *u« Pom 1300*34 unidiretfiona, lamma prope , 

at -250°F usinn the CCA formulation of Hashin lor predicting effective elastic 

: " , , „„ r.TE data using the Rosen and Hashin analysis 

and +250 °F from T300/934 unidirectional lamina C 

discussed in Appendix A. 
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Table 2. Matrix Properties at 75 °F . 


Matrix 

E 

(Msi) 

g' 1 ) 

(Msi) 

1/ 

a 

(10-®/ 0 F) 

934 epoxy< 2 > 

0.67 

0.25 

0.36 

23.0 

5208 epoxy' 2 * 

0.67 

0.25 

0.36 

23 0 

930 epoxy' 2 ) 

0.67 

0.25 

0.36 

23.0 

CE339 epoxy' 2 ) 

0.67 

0.25 

0.36 

35.2< 3 > 

PMR15 polyimide 

0.50 

0.19 

0.35 

20.0 

2024 aluminum 

10.60 

4.00 

0.33 

12.9 

Borosilicate glass 

9.10 

3.80 

0.20 

1.8 


(1) G = E/2(1 + u) 


(2) All epoxies were assumed to have the same properties unless noted 

(3) Value for a similar rubber-toughened epoxy 


Second order polynomial representations of the dependence of the 934 properties with tempera- 
ture are given in Table 3. These polynomials were used in the analyses for computing the matrix 
properties at discrete temperatures between -150°F and +250°F 


Table 3. Functional Form of Temperature Dependent 934 Epoxy Properties. 


Property = A 0 + A 1 T + A 2 T 2 , (T = °F) 


Property 

Ao 

Ai 

A2 

E (Msi) 

0.7743x10° 

-0. 1422X10- 2 

0.1123x10 -5 

G (Msi) 

0.2832x10° 

-0.51 60x1 0~ 3 

0.4448x10-° 

i/ 

0.3677x10° 

-0.4400x1 0“ 4 

-0.2514x10-° 

a (10 _6 /°F ) 

0.1 975x1 0 2 

0.4500x1 O’ 1 

-0.1 700x1 0" 4 


5.2 Effective Coefficients of Thermal Expansion 


Predictions of effective CTE's from the four explicit formulations presented in Section 2.1, and 
from the finite element analysis, developed for this study and presented in Section 3.1.1, will be 
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compared with each other and with experimental data. These five analyses will be referred to by 
the following abbreviations: Shapery (SH), Chamberlain (CB), Chamis (CH), Rosen and Hashin 
(RH), and finite element (FE). The Hex and Sq designations for the CB and FE predictions refer to 
hexagonal and square fiber arrays. 

The specific material systems examined are given in Table 4, along with their experimentally 
determined values of a 1 and a 2 at RT. Some of these values were measured specifically for this 
study, while others were taken from previous studies* 61 -62) . All of the data were obtained from the 
interferometric dilatometer system described in Section 4.1. These material systems have axial 
fiber to matrix stiffness ratios of E f t /E m ranging from 6 to 140, and axial fiber to matrix CTE ratios of 
a\ /a m ranging from -0.01 to -0.30, and thus cover a wide range of fiber/matrix combinations. The 
RT fiber and matrix properties needed as input to the analyses were given in Tables 2 and 3. 

Predictions and comparisons with experimental data for the axial CTE will be discussed first, 
followed by a discussion of transverse CTE predictions and comparisons. These predictions will 
assume temperature independent constituent property behavior. The sensitivity of effective CTE’s 
to constituent properties, the effects of temperature dependent constituent properties, and the 
effects of fiber orthotropy will be discussed in seperate sections. 


Table 4. Experimentally Determined CTE Values at 75 °F . 


Material System 
(Fiber/Matrix) 

V* 

(%) 

«1 

(10- 6 /°F) 

O 

O 

T300/5208 

68 

-0.063 

14.02 

T300/934 

57 

-0.001 

16.13 

P75/934 

48 

-0.584 

19.18 

P 75/930 

65 

-0.598 

17.62 

P75/CE339 

54 

-0.567 

26.34 

C6000/PMR15 

63 

-0.118 

12.46 

HMS/Borosilicate 

47 

-0.230 

2.10 

PI 00/2024 

40 

0.800 

14.51 
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5.2.1 Axial Coefficient of Thermal Expansion 

Figures 11-14 show a comparison of a i as a function of V 1 predicted from the different analyses for 
four of the material systems: T300/934, P75/934, HMS/Borosilcate, and PI 00/2024. Experimental 
data are also shown on these figures. Predicted and measured values of ai for all of the materials 
studied are given in Table 5. The SH, CB, and CH analyses all used equation (2.1 .1) for predicting 
ai and are labeled SH in the figures and Table. As shown in the figures, the differences between 
the RH and FE results were negligible and there were only small differences between these and 
the SH results. All of the analyses were in good agreement with the experimental data The largest 
deviation between any of the predicted and experimental values was only 0.1 x 10~ 6 /°F , and in 
most cases the deviation was on the order of about 0.05x 10- 6 /°F • Although the magnitudes of 
differed for the different material systems, the general response was the same (i.e. decreasing 
a, with increasing V*). This implies that the relative magnitudes of the fiber/matrix stiffness and 
CTE ratios did not significantly affect the general trend in at as a function of V*. This was not true 
for 02 , as will be discussed subsequently. 

Table 5. Comparison of Experimental and Predicted Values of the Axial CTE at 75 F . 




«,(io-® 



Material 





FE 

System 

Exp 

SH 

RH 

Hex 

Sq 

T300/5208 

-0.063 

-0.083 

-0.051 

-0.050 

-0.041 

T300/934 

-0.001 

0.045 

0.089 

0.089 

0.093 

P75/934 

-0.584 

-0.535 

-0.511 

-0.512 

-0.511 

P75/930 

-0.598 

-0.642 

-0.627 

-0.627 

-0.623 

P75/CE339 

-0.567 

-0.493 

-0.460 

-0.461 

-0.460 

C6000/PMR1 5 

-0.118 

-0.125 

-0.104 

-0.104 

-0.099 

HMS/Borosilicate 

-0.230 

-0.180 

-0.180 

-0.181 

-0.181 

PI 00/2024 

0.800 

0.875 

0.907 

0.906 

0.905 
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Analysis Experimental data 


SH.CB. CH 
RH, FE 


• T300/934 
■ T 300/ 5208 



Figure 11. Axial CTE of unldlr.c.lon.1 T30W934 end 7300/5208. 


Analysis Experimental data 



Figure 12. Axial CTE ol unidirectional P75/934 and P75/930. 
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• Experimental data 



Analysis 

SH, CB, CH 
RH, FE 


Figure 13. Axial CTE of unidirectional HMS/Borosllicate. 



Analysis 

SH, CB, CH 
RH, FE 


Figure 14. Axial CTE of unidirectional PI 00/2024. 
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5.2.2 Transverse Coefficient of Thermal Expansion 


Figures 15-18 show the response of a 2 as a function of V* for the same four material systems shown 
in figures 10-13. Experimental data are also shown on these figures, as well as in Table 6, where 
predicted and experimental values are compared for all of the material systems. The SH results 
were computed with the modified version for transversely isotropic fibers (equation (2.1.3)). 

There were large differences between the predicted values of all of the analyses, except for the 
RH and FE analyses. Results from these two analyses were in excellent agreement with each 
other, including both array geometries of the FE model. Significant differences between the CH 
and FE analyses for prediciting « 2 have been previously documented in the literature* 20 *. These 
differences were attributed to Poisson restraining effects not included in equation (2.1.6). The 
neglection of this type of three dimensional effect was also thought to be responsible for the large 
difference between the CB and FE results. Figures 15-18 and Table 6 also demonstrate that the RH 
and FE results for n 2 were consistently in much better agreement with the experimental data than 
were the other analyses. The SH results were in slightly better agreement with the experimental data 
for the P75/930 and P75/CE339 systems. However, it should be noted that the matrix mechanical 
properties for these two systems were assumed to be the same as the other epoxy matrices, which 
is probably not an accurate assumption, and therefore the better agreement with the experimental 
data is believed to be fortuitous. It should also be remembered that the modification of the SH 
analysis for transversely isotropic fibers was not based on any mathematical derivation, and was 
included for comparison purposes only. Agreement between experimental values and RH and FE 
„oicted values were usually within about 1 5 percent. Predictions from the other analyses differed 
with the experimental data by as much as 50 percent. 
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Experimental data 


Analysi 



0 20 40 - 60 80 100 

V f (%) 


Figure 15. Transverse CTE of unidirectional T300/934 and T30/5208. 


Experimental data Analysis 



Figure 16. Transverse CTE of unidirectional P75/934 and P75/930. 
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RH 



Figure 17. Transverse CTE of unidirectional HMS/Boroslllcate. 


• Experimental data Analysis 



Figure 18. Transverse CTE of unidirectional PI 00/2024. 
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Table 6. Comparison of Experimental and Predicted Values of the Transverse CTE at 75 °F 


Material 

System 




a 2 (10- 

- 6 /°F) 




Exp 

SH 

CB 

CH 

RH 

FE 

Hex 

Sq 

Hex 

Sq 

T300/5208 

14.02 

14.60 

8.85 

7.28 

10.10 

13.00 

12.99 

12.72 

T300/934 

16.13 

17.30 

10.60 

9.20 

11.90 

15.60 

15.69 

15.59 

P75/934 

19.18 

18.60 

11.10 

9.63 

12.20 

17.90 

18.05 

18.03 

P75/930 

17.62 

14.10 

7.89 

6.20 

9.13 

13.20 

13.37 

13.33 

P75/CE339 

26.34 

24.70 

13.80 

11.30 

15.50 

23.60 

23.94 

23.88 

C6000/PMR15 

12.46 

14.30 

8.92 

7.67 

10.00 

1240 

12.56 

12.63 

HMS/Borosilicate 

2.10 

3.32 

3.04 

3.18 

3.27 

2.49 

2.46 

2.47 

PI 00/2024 

14.51 

11.90 

8.08 

7.47 

8.52 

15.00 

15.15 

15.13 


Unlike the results for aj , the response of a 2 as a function of V* was affected by the fiber/matrix 
property ratios. T300/934 and P75/934 had similar fiber/matrix property ratios, and exhibited a 
similar response (Figures 1 5-1 6). The PI 00/2024 system (Figure 1 8) had CTE fiber/matrix property 
ratios similar to T300/934 and P75/934, but had much smaller moduli ratios. This difference in 
moduli ratios resulted in a different response of a 2 as a function of V* as predicted from the RH and 
FE analyses (i.e. increasing a 2 with increasing V 1 for small values of V f ). The HMS/Borosilicate 
system (Figure 17) had moduli ratios similar to PI 00/2024 Al but had much larger CTE property 
ratios. This difference caused a significantly different response in a 2 as a function of V* compared 

to the trends exhibited in Figures 15-16 and 18 (i.e. increasing a 2 with increasing V* for all values 
of V*). 


The results presented above and from Section 5.2.1 show that both simple and rigorous analyses 
accurately predict the axial CTE of carbon fiber unidirectional composites. These results also 
demonstrate that plane stress (CB analysis) and strength of materials (CH analysis) formulations 
do not accurately predict the transverse CTE of these same composites. The more rigorous RH 
analysis does accurately predict the transverse CTE and can be used in lieu of a detailed FE 
analysis for such predictions. 
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5.2.3 Sensitivity to Constituent Properties 


As stated earlier, the analytical predictions discussed in the previous section used constituent 
property data compiled from various literature sources. Some of the values were measured directly, 
some were back-calculated from composite data, and some were only estimates. A parametric 
study was conducted using the RH analysis to determine the sensitivity of o, and a 2 of P75/934 
to the constituent properties. Each of the independent thermoelastic fiber and matrix properties 
were individually reduced 20 percent, while holding the other properties constant. The resulting 
changes in «i and c* 2 are 9' ven ' n Table 7. 

As would be expected, ay was most sensitive to E f , and . Fortunately, these are the easiest fiber 
properties to measure, and therefore usually the most reliable. The prediction of ay was insensitive 
to the transverse fiber properties and matrix properties. The prediction of a 2 was most sensitive 
to matrix properties. Matrix properties are relatively easy to obtain experimentally, however there 
are concerns as to whether the chemistry of resin matrices is the same in bulk form as it is in the 
composite. The sensitivity of a 2 to transverse fiber properties was less than expected. This is again 
fortunate because transverse fiber properties are the most difficult to obtain, and therefore usually 
the most unreliable. 

The sensitivity analysis presented above demonstrates that the good agreement between the RH 
and FE analyses and the experimental data was not fortuitous. The constituent properties with the 
most uncertainty (transverse fiber properties) have only a small effect on the prediction of a, and 

a 2 . 

The sensitivities of a! and a 2 to V* for all of the material systems studied have been illustrated in 
Figures 11-18. Predicted values of at were very sensitive to small changes in V* for values of V 1 
less than 30 percent. This sensitivity decreased for large values of V* . The sensitivity of a 2 to V* is 
much less, from a percentage standpoint, than that of 
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Table 7. Sensitivity Analysis for P75/934 CTE’s (V f =0.60). 


20 Percent Reduction in 

Percent Change 

a 2 

| 

-6.2 

+0.2 

-0.6 

+5.8 

g ; 2 

gL 

0.0 

0.0 

0.0 

-5.2 

4 

i 

E* 

-0.7 

0.0 

-24.9 

-0.6 

0.0 

-3.2 

+ 5.0 

-10.4 


+ 1.5 

-7.1 

a m 

+5.0 

-16.9 


(1) Calculations based on Rosen and Hashin Analysis 


5.2.4 Effects of Temperature Dependent Constituent Properties 

As discussed in Section 5.1, the 934 epoxy matrix properties are not temperature independent. 
The method for properly incorporating the effects of temperature dependent constituent properties 
in the prediction of effective CTE’s was described in Section 3.2. This method was applied to the 
RH analysis, and the results will be presented in this section. The RH analysis was selected due to 
its consistently better agreement with experimental data than the other closed-form analyses, and 
because of its computational simplicity compared to the FE analysis. 

Figures 19-21 show temperature dependent CTE predictions, labeled exact, compared with ex- 
perimental data for T300/934, T300/5208, and P75/934. Also shown in each figure are predictions 
from the RH analysis using an approximate method for incorporating temperature dependent matrix 
properties (i.e. using a m at temperature T, rather than the a m dT). The exact analysis predic- 
tions were in better agreement with the experimental values of a 2 than the approximate analysis 
for all three material systems. The excellent agreement between predicted values and experimen- 
tal data for a 2 of T300/934 (Figure 1 9) were expected, due to the fact that the T300/934 data was 
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CTE (° 


Figure 19. 


CTE ( 


Figure 20. 


E xac t l Analysis 
Approximate! 



Temperature dependent CTE of unidirectional T300/934 (V*— 0.57). 



Temperature dependent CTE of unidirectional T300/5208 (V*— 0.68). 


Results and Discussion 


47 




Figure 21. Temperature dependent CTE of unidirectional P75/934 (V^=0.48). 

used to back-calculate the temperature dependent properties of the 934 matrix. These figures do 
demonstrate the error introduced with the approximate analysis. The agreement between predicted 
and experimental temperature dependent values of was also good. 

The above results demonstrate that the CTE’s of carbon fiber reinforced epoxies do vary with 
temperature, and accurate predictions of these values require accurate data on the temperature 
dependence of the constituent properties. 


5.2.5 Effects of Fiber Orthotropy 

As mentioned previously, Knott* 52 ) conducted an extensive analytical study of the effects of carbon 
fiber orthotropy on the effective properties and stresses in unidirectional composites, using the 
CC model discussed in Section 3.3.2. The effective CTE's predicted with the assumption of 
transversely isotropic (Tl) fibers were compared to those obtained by assuming two types of 
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cylindrical fiber orthotropy; circumferential, where Eg > E r , and radial, where E r > Eg. For 
transversely isotropic fibers, Eg = E r . Large differences in the effective CTE’s were predicted 
for the three types of fiber orthotropy. However, the fiber properties used in that study represent 
extreme degrees of cylindrical orthotropy. For the circumferentially orthotropic (CO) case E fl and 
a e were set equal to E* and o x , respectively, while E f and a r were maintained at their transversely 
isotropic values. For the radially orthotropic (RO) case E r and a, were set equal to E x and a x , 
respectively, while Eg and Q0 were maintained at their transversely isotropic values. No data in the 
literature could be found to verify that this degree of cylindrical orthotropy actually exists in carbon 
fibers. As mentioned earlier, and as pointed out by Knott, the experimental data suggests that 
carbon fibers contain a transversely isotropic core, and therefore the effective degree of cylindrical 
orthotropy for the entire fiber would be much less than values described above. 

The effects of a less severe degree of cylindrical orthotropy were considered for the present study. 
The fiber properties for T300 and P75, given earlier in Table 1, were modified according to the 
relations shown in Table 8 to account for cylindrical orthotropy. The CC model was used to 
compute the effective CTE’s for the three types of fiber orthotropy shown in Table 8. As described 
in Section 3.3.2, this model is formulated on the basis of cylindrical material orthotropy, where 
transverse isotropy is a special case of the more general condition. The results for the CTE s of 
T300/934 (V f = .57) are shown in Figure 22, and the results for the CTE s of P75/934 (V* = .48) 
are shown in Figure 23. Also shown on these figures are predictions from the RH analysis for 
transversely isotropic fibers. There are two main points demonstrated by these figures. First, the 
RH analysis and the CC model predict very similar values for the case of transversely isotropic 
fibers. This implies that the CC model is another formulation, in addtion to the RH and FE 
discussed in Sections 5.2.1 and 5.2.2, that can be used to accurately predict the effective CTE’s 
of unidirectional composites. Secondly, the degree of cylindrical orthotropy shown in Table 8 


Results and Discussion 


49 



produces only small differences, in absolute magnitudes, in predicted CTE’s compared to the 
assumption of transversely isotropic fibers. 


Table 8. Cyllndrically Orthotropic Carbon Fiber Properties. 


Property 

Circumferential 

Orthotropy 

Radial 

Orthotropy 

Transverse (,) 

Isotropy 

E f 

= 2E| 

= 2E** 


G* 

G'x* = 2G* r 

G!cr = 2G f xfl 

Gjcr = G f x , 


v \.9 = *®*4r 

^xr = 

*4 = "le 

a* 

Q e ~ • 5a< r 

Q r = 

q\ = a\ 


(1) x, r, and 9 directions are equivalent to the 1, 2, and 3 directions, respectively. 


Based on the above results, carbon fiber transverse isotropy is a reasonable and simplfying 
assumption for predicting effective CTE’s of unidirectional carbon fiber reinforced composites. This 
is especially true given the lack of quantitative experimental data on fiber mechanical and thermal 
properties. 

5.3 Thermally Induced Stresses 

The majority of thermal stress analysis results are presented for P75/934. This is a high modulus 
composite system with many potential applications on stiffness and dimensionally critical space- 
craft structures. These results were generated with the composite cylinder (CC) and finite element 
(FE) models described in Section 3.3, and used the fiber properties from Table 1 and the tem- 
perature dependent 934 matrix properties from Table 3. In the analyses, the principal material 
coordinates (1,2,3) of Figure 4 correspond to the FE cartesian (x,y,z) coordinate system of Figure 
5, and the CC cylindrical (x,r,0) coordinate system of Figure 7. Stress components are presented in 


Results and Discussion 


50 



□ Tl, RH analysis 
CO > . CC analysis 



«1 “2 

Figure 22. Effects of fiber orthotropy on the CTE of T300/934 (^=0.57). 


IZD Tl, RH analysis 



«l «2 

Figure 23. Effects of fiber orthotropy on the CTE of P75/934 (^=0.48). 
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terms of the cylindrical coordinate system (i.e. <r x , <r r , <rg , and r r0 ). The thermal load used, unless 
otherwise noted, was a uniform AT of -500 °F . This AT corresponds to a stress- free-temperature 
(SFT) of +350 °F , which is the cure temperature for 934 epoxy, and a use temperature of -150 F , 
a typical spacecraft operating environment. Results for unidirectional laminates will be presented 
first, followed by results from the global/local formulation for predicting stresses in multidirectional 

laminates. 

5.3.1 Unidirectional Laminates 

5.3.1. 1 Comparison of Analyses. Thermal stress results from the CC and FE analyses 
are presented for unidirectional P75/934, with a fiber volume fraction, V*. of 0.60, in Figures 24-29. 
The matrix stresses at the fiber/matrix interface as a function of the circumferential position around 
the fiber ( 0 is measured from the y axis of Figure 5) are shown in Figures 24 and 25. The maximum 
values of <t x (16.2 Ksi) and a e (15.2 Ksi) occurred on this interface at the circumferential locations 
shown in the figures. Both the CC and FE analyses predict approximately the same magnitude and 
distribution for a* . The FE results for <r 9 do exhibit a small dependence on 0, and differ in magnitude 
from the CC results. For all cases, the FE results using the hexagonal fiber array geometry were in 
closer agreement with the CC results than those using the square array. This was expected, due 
to the fact that the hexagonal array more closely resembles an axisymmetric geometry. 

The largest differences between the analyses are exhibited in the <r r and t , 0 components at the 
interface. As shown in Figure 24, the a r component varied from approximately -7 Ksi at 0 = 0°, 
to a maximum tensile value of 0.60 Ksi at approximately 0 = 37° for the square array FE results. 
The CC results predicted a constant value of -3.3 Ksi. There was also a significant r re component 
(Figure 25) in the square array FE results that was nonexistent in the CC results. The maximum 
value of r, e was 5 Ksi (not shown on Figure 25), which occurred at 0 = 22.5°, not on the interface 
but at r/r f = 1 .2, where r f is the fiber radius. 
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[0] P75/934, V f = 0.6, AT = -500° F 



Figure 24. Axial and radial Interfaclal matrix stresses. 


[0] P75/934, V f =0.6, AT = -500°F 



Figure 25. Tangential and shear Interfaclal matrix stresses. 
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[0] P75/934, V* = 0.6, AT = -500°F 



Normalized radial coordinate, r/r f 


Figure 26. Tangential stresses along the line of maximum fiber spacing. 


[0] P75/934, V f = 0.6, AT = -500° F 



Figure 27. Axial and radial stresses along the line of maximum fiber spacing. 
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[0] P75/934, V f = 0.6, AT = -500° F 



Figure 28. Tangential stresses along the line of minimum fiber spacing. 


[0] P75/934, V f = 0.6, AT = -500° F 



Figure 29. Axial and radial stresses along the line of minimum fiber spacing. 
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Figures 26 and 27 show the stresses along a radial line of maximum distance between fibers. The 
circumferential position of this line is either 9 = 30°, or 45° depending on whether the results are 
for the hexagonal or square array, respectively. The radial coordinate was normalized with respect 
to the fiber radius. The CC results are the same for all radial lines due to its axisymmetric solution. 
There is no t , 0 predicted along this radial line with the FE analysis because it is a line of symmetry. 
Both the CC and FE analyses predict similar magnitudes and distributions for the fiber stresses. All 
three components are compressive, with a maximum value of approximately -1 1 Ksi for <r x . Both 
analyses also predict approximately the same magnitude and ditribution for the a % matrix stress. 

There are significant differences between the predicted values of a, and a 0 in the matrix from the 
two analyses. The a 0 component predicted from the FE analysis is significantly smaller than <jq 
from the CC analysis (Figure 26). The other major difference between the two analyses is exhibited 
in the <r r component along this radial line (Figure 27). The CC analysis predicts a <r r stress that 
is compressive at the interface and decreases to zero at the outer boundary of the CC model, as 
required by the free surface boundary conditions. The FE analysis also predicts a compressive a, 
at the interface, but this component reaches a tensile value of 5.3 Ksi at the boundary of the FE 
model (i.e the midpoint of maximum fiber spacing). This is because the FE model does not treat 
the fiber and surrounding matrix as an isolated problem with stress free boundaries, but rather 
takes into account the influence of adjacent fibers. 

Figures 28 and 29 show the stresses along the radial line of minimum distance between fibers 
(i.e. 0 = 0°). The trends are similar to those exhibited in Figures 26 and 27. The CC results are 
identical to those shown in Figures 26 and 27 due to the axisymmetric nature of its solution. Again, 
the shear stress along this line is zero, and the fiber stresses are compressive. The most significant 
differences between the stresses along this line and those along the line of maximum fiber spacing 
occur in the matrix n g and a, components from the FE analysis. The interfacial value of the og 
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component (Figure 28) is significantly larger at this location (15 Ksi) compared to the interlacial 
value at the maximum distance between fibers (1 1 Ksi. Figure 26). Also, the <r r component (Figure 
29) remains compressive in the matrix with a magnitude similar to that of the fiber <r r component. 

The above results show that FE analysis predicts higher localized matrix stresses than the CC 
analysis. This is because the FE analysis takes into account the influence of adjacent fibers, rather 
than treating a single isolated fiber as in the CC analysis. The effects of these localized higher 
stress regions on the initiation of damage will be discussed in a subsequent section. 

5.3.1 .2 Effects Of Constituent Properties. A comparison of the stress distributions for 
P75/934 and T300/934 are shown in Figure 30. The normal stresses at the interface are compared 
for the FE square array analysis at a V» = .60. The differences in the stresses are relatively small 
compared to the large differences in fiber moduli and CTE (Table 1). This may be explained by 
examining a simple one-dimensional strength of materials formulation for predicting axial thermally 
induced stresses. Assuming temperature independent properties, the axial stress in the fiber and 

matrix may be written as (63> 


4 = v m E m E 


, ( 

1 yv*E , 1 + v m 


V nr 'E m 


AT 


m _ _v f E m E f ( — — — — ^ /IT 

~ V E ^v f E* + V m E m ) 


(5 3.1.2.1) 


(5.3.1. 2.2) 


where V m = (1 -V f ). In polymer matrix composites reinforced with graphite fibers, where E^ »E , 
\Z<e* + V m E m may be approximated by V*E!| . This approximation results in less than a 2% error 
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[0] FE-Sq analysis, V f = 0.6, AT = -500° F 



0, deg 

Figure 30. Effects of fiber modulus on Interfacial matrix stresses. 

fork's greater than .50. The term (« m - a\) may also be approximated by o m without introducing 
greater than a 3% error. This is because r* m » rv^ . 

These two approximations result in simplified forms of equations (5.3.1. 2.1) and (5.3.1 .2.2) written 


<r f x = (V m /V , )E m rv"MT (5.3.1. 2.3) 

and 

< = -E m « m /lT (5.3.1 .2.4) 

The fiber properties do not appear at all in equations (5.3. 1.2. 3) and (5.3. 1.2.4). Although the above 
simplified analysis is based only on the axial stress component, the same trend is exhibited by the 
other components, as evidenced by Figure 30. It is also interesting to note that V* does not appear 
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in equation (5.3.1 .2.4), implying that the axial matrix stresses are not a strong function of V f for V f, s 
greater than .50. 

In order to significantly reduce the magnitudes of thermal stresses for a fixed AT , the term E m a m 
must be reduced. Unfortunately, there exsists an inverse proportionality between E m and a m for 
a wide range of polymer systems^ 64 ), as shown in Figure 31. The product, E m a m , ranges from a 
high of 1 7.7 to a low of 1 2 psi/°F for the polymers shown in the figure. The 934 epoxy has a value 
of 15.4 psi/°F . Therefore, for a given AT , the matrix stresses in a 934 system could be reduced 
by approximately 22% by using one of the other polymers shown in the figure (i.e. polyphenylene) 
as the matrix. However, it is not known whether all of these other polymers would make suitable 
matrix materials for graphite reinforced composites. Large reductions in the level of residual thermal 
stresses in graphite reinforced composites are obtained by lowering the applied AT during use. 
This could be accomplished from a materials standpoint by lowering the SFT through processing 
and/or chemistry modifications. 

The importance of matrix properties on thermally induced stresses was demonstrated above. 
Figure 32 shows the importance of properly accounting for the temperature dependence of these 
properties. Matrix stresses at the fiber/matrix interface were compared for the cases of temperature 
independent matrix properties and temperature dependent matrix properties. The RT values of 
934 epoxy shown in Table 2 were used for the temperature independent case. These results 
were generated with the FE analysis assuming a square array of fibers. As shown in the figure, 
the temperature independent property results exhibit stresses that are significantly smaller than 
the te*~ mature dependent property results for a x and ctq by as much as 38%. This is due to 
the stiffening of the 934 epoxy matrix at low temperatures which is neglected when temperature 
independent properties are assumed. 

Based on the results presented above, the effects of fiber orthotropy on thermally induced stresses 
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O Epoxy(934) 

□ Polystyrene 
O Polyester 
Polycarbonate 
A Lexan 101 
A Lexan 3414 
Polyphenylene 
(2i Polyimide 

Polyethersulfone 
Q> Victrex 200T/300T 

* Victrex 420T/520T 
Q PEEK 

Polyamide-imide 
O Torlon 4203 

♦ Torlon 9040 
0 ABS/PVC 


Figure 31. Relationship between CTE and modulus for various polymers. 


[0] P75/934, FE-Sq analysis, V f = 0.6, AT = -500°F 



Figure 32. Effects of temperature dependent matrix properties on interfacial matrix stresses. 
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would be expected to be small. This was proven to be true, as evidenced by the results ptesented 
in Table 9. These results were generated withthe CC analysis, and used the cylindrical^ orthotropic 
properties shown in Table 8. Only small ditferences in the maximum matrix stresses were predicted 
lor the three types of liber orthotropy; transversely isotropic fn). radially orthotropic (BO), and 
circumferentially orthotropic (CO), The RO fiber had the largest predicted deviation from the Tl 
case, with an approximately 11% larger a e stress component. 

Table 9. Effects of Fiber Orthotropy on P7S/934 Thermal Stresses <V=0.60, /ST=-500"F ). 


Maximum Matrix Stress* 11 (Ksi) 


Type of Orthotropy 


< U 

*6 

Tl 

15.96 

-3.33 

13.31 

RO 

16.36 

-3.70 

14.80 

CO 

15.99 

-3.36 

13.43 


(1) Calculations based on CC Analysis 


5.3.1 .3 Effects of Fiber Volume Fraction. The effects of V on the stress distributions in 
unidirectional P75/934 were also determined for each analysis. The rr x , rr, , and tr e matrix stress 
components, at the fiber/matrix interface, are shown as a function of V in Figures 33, 34, and 
35, respectively. These results were generated with the FE analysis assuming a square array of 
fibers. The results for the hexagonal array show similar trends. In general, the FE analysis pre- 
dicted increases in the absolute value of all of the stress components with increasing V. The only 
exception to this was the *, component at 0 = 45°. which changed from a compressive -5.7 Ksi 
a V - . 1 0 to a tensile value of 2.8 Ksi at V' = .75. The behavior of the matrix c, component was 
also the major u„lerence between the stresses as a function V< from the CC and FE analyses The 
analysis predicted a steadily decreasing compressive value of m towards zero with increasing 
V. One of the most significant features of these daa is the presence of large matrix stresses even 
at very low values of V*. This is especially true for rr x , where the maximum matrix stress at V' - 10 
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[0] P75/934, FE-Sq analysis, AT = -500°F 
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Figure 33. Effects of V* on axial interfacial matrix stresses. 


[0] P75/934, FE-Sq analysis, AT = -500°F 
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Figure 34. Effects of V* on radial Interfacial matrix stresses. 
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[0] P75/934, FE-Sq analysis, AT = -500° F 
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Figure 35. Effects of V* on tangential interfacial matrix stresses. 


is only 30% less than the maximum value of 17.6 Ksi at V* = .75. These results support the trends 
exhibited by the simple one-dimensional analysis presented in the preceding section, where V* 
does not appear in the expression for (equation (5.3. 1.2.4)). 


The array geometries used in the FE analysis have upper bounds on V* for contiguous fibers. These 
upper bounds are V f = .785 and V* = .907 for the square and hexagonal arrays, respectively. 
The upper bound on V* for the CC analysis is 1.0, or no matrix phase, for which case the matrix 
stresses are obviously zero. The matrix stresses at the fiber/matrix interface for V f ’s of .75, .85, and 
.99, corresponding to the FE square array, FE hexagonal array, and CC analyses, respectively, are 
shown in Figure 36 and 37. Values of .75 and .85 were used as upper limits on V* for the square 
,d hexagonal array FE analyses, respectively, due to difficulties in modeling a contiguous fiber 
array with the type of elements used in the analysis. An upper limit of V* = .99 was used for the 
CC analysis to represent the presence of a small, but measurable matrix phase. It is interesting to 
note that, for the a x and <r e components, there are only small differences in the maximum tensile 
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[0] P75/934, AT = -500 °F 



Figure 36. 


Effects of analysis V* upper bound on 


axial and radial Interfaclal matrix stresses. 


[0] P75/934, AT = -500°F 



Figure 37. 


Effects of .navel. V upper bound on tangential end .hear lnf«rf«l.l matrix 
stresses. 


Results and Discussion 


values for the three different values of V f . 

A real composite does not have fibers arranged in regular periodic arrays, as assumed in the FE 
analysis, nor consists of cylinders of fiber and matrix with varying diameters, as assumed in the 
CC analysis. A photomicrograph of the structure of a real composite is shown in Figure 38. One of 
the most striking features revealed by this figure is the wide range of V* ’s present in the specimen. 
Even though this specimen has an average or global V* of .50, there are localized regions with 
much smaller and much larger N^’s. 

An attempt at modeling localized regions of high V* in a global V* region of 0.60 was made using 
a global/local FE scheme, similar to the one described in Section 3.3.3. The procedure is shown 
schematically in Figure 39. The FE model for V< = .75 in the square array was constrained by 
three different sets of global boundary conditions as follows: Case 1, < x and e y from V< = 60, 
Case 2, c x from V f = .60 and f y unconstrained, and Case 3, e x from V f = .60 and e y from V* = .75. 
The results from this procedure, for the matrix stress components at the fiber/matrix interface, are 
shown in Figures 40 and 41 . Cases 2 and 3 gave nearly identical results and therefore are labled 
as one curve in the figures. For comparison purposes, the stresses for actual cases of V* = .60 
and V* = .75 are also shown in the figures. 

As can be seen in the figures, cases 2 and 3 gave results very similar to the actual V* = .75 case. 
Case 1 predicted stresses lower than those of the actual V 1 = .60 case. It should be noted that 
cases 1 , 2, and 3 all predicted tensile a x fiber stresses. This is opposite of the compressive <r x fiber 
stresses predicted in all of the earlier results (Figures 26-29), and is due to the imposition of the 
V = .60 fx boundary condition. The V f = .75 model normally contracts more in the axial direction 
than is allowed by this boundary condition, and thus tensile axial stresses are set up in both the fiber 
and matrix. Although none of the three cases accurately represent the real boundary conditions, 
cases 2 and 3 seem more plausible since they maintain a uniform axial strain consistent with the 
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Figure 38. Photomicrograph of typical composite showing local variations In V*. 




e y (V» = .60) 


e x (Vf = .60) 
Case 1 



£y (V* = .75) 


Case 2 


Case 3 


Figure 39. Global/local finite element analysis of localized high V* regions. 


Results and Discussion 


[0] P75/934, FE-Sq analysis, AT = -500°F , 



0, deg 


Figure 40. Effects of high local V* on axial and radial Interfaclal matrix stresses. 


[0] P75/934, FE-Sq analysis, AT = -500° F 



0, deg 


Figure 41. Effects of high local V 1 on tangential and shear Interfaclal matrix stresses. 
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global V f , but allow the transverse strains to be consistent with the local V f . Case 1 seems overly 
constrained by forcing the transverse strains to be consistent with the global V* over a distance of 
one fiber radius. The real boundary conditions probably lie somewhere between cases 2 and 3, and 
case 1 . Based on these results, it appears that a resonable conservative first approximation to the 
matrix stresses in localized regions of high V* can be predicted from the actual (i.e unconstrained) 
model of that region. 

5.3.2 Multidirectional Laminates 

The global/local formulation, described in Section 3.3.3, was used to determine the thermally 
induced matrix stresses in multidirectional P75/934 laminates. The class of laminates selected for 
investigation were the [0 2 / ±0] s family, with 0 having values of 0, 1 5, 30, 45, 60, 75, and 90 degrees. 
This class of laminates includes both the minimum amount of ply property mismatch (0 = 0°, or 
unidirectional), and the maximum amount of ply property mismatch ( 6 = 90°), or “worst case”. 
Lamina strains in the principal material coordinate directions were determined for these laminates 
using equations (3.3.3. 1) through (3.3.3.5). The lamina elastic properties used in these calculations 
were determined from the CCA model of Hashin< 23 > using the fiber and matrix properties given in 
Tables 1 and 3. The calculated lamina strains for a AT of -500 °F are shown in Table 10. These 
values were used as boundary conditions for the FE micromechanics analysis. 

A comparison of unidirectional and [O 2 / ± 0] s matrix stresses in the ±9 plies at the fiber/matrix 
interface are shown in Figures 42-44. The <r x , . <?o < and r r g components (Figures 42-43) all 

increase in magnitude with increasing ply angle up to an angle of ±60°. For angles greater than 
±60° the magnitudes of the stresses remained nearly constant. These layups are labeled as one 
curve in the figures. The most significant difference between the unidirectional and [O 2 / ± 0]s 
response is the presence of a large tensile value of a , at the fiber/matrix interface. The a, 
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component reaches a maximum interfacial tensile value of 7 Ksi in the [ 02 / 902 ]* laminate at a 
circumferential position of approximately 30°. The overall maximum tensile value of the matrix <r r 
component occurs at this circumferential position, but at the midpoint between fibers. The other 
major difference between the unidirectional and [O 2 / ± 0]s response is the presence of out-of- 
plane shear stresses t#* and t X i*. The maximum matrix values of these components occurred at 
the fiber/matrix interface, and are shown in Figure 44. Except for the presence of these out-of- 
plane shear components, the magnitudes and trends of the matrix stresses in the 0° plies were 
very similar. 


Table 10. Lamina Strains in [0 2 /±0]* P75/934 Laminates (^=0.60, ziT=-500°F ). 


Laminate 

Orientation 

0° 

Ply Strains* 1 ) 

no- 6 ) 

±6 

Ply Strains* 1 > 

(io- 6 ) 

<1 

e 2 

712 

*3 

<4 


712 


[0] 

261 

-7673 

0 

-7800 

na 

na 

na 

na 

[Os/ ±15]s 

488 

-7180 

0 

-8200 

-26 

-6670 

-3830 

-8200 

[0 2 /±30] 8 

654 

-3730 

0 

-10000 

-442 

-2640 

-3800 

-10000 

[Os/ ± 45]* 

358 

-991 

0 

-11100 

-316 

-316 

-1350 

-11100 

[0 2 / ± 60]* 

132 

-214 

0 

-11400 

-128 

45 

-300 

-11400 

[0 2 / ± 75]* 

28 

-34 

0 

-11400 

-30 

24 

-31 

-11400 

[0 2 /90 2 ] s 

~ 0 

- 0 

0 

-11400 

~0 

~ 0 

0 

-11400 


(1) Calculations based on LPT 


There were also differences in the fiber stresses between the [O 2 /±0 ]s and unidirectional laminates. 
The largest difference was observed in the <r x component. The [0 2 / ± 15] s , [0 2 / ± 30] s , and 
[0 2 / ± 45] s laminates all had predicted tensile values of the fiber <r x component in the 0° plies. The 
v.yest of these was 74 Ksi in the [0 2 / ± 30] s laminate. This was in contrast to the compressive 
value of approximately -1 0 Ksi predicted for the unidirectional laminate. The other major difference 
was the -100 Ksi <r x fiber stress in the ±15° plies of the [0^ ± 15]* laminate. This value was an 
order of magnitude larger than the value predicted in the unidirectional laminate. 
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P75/934, FE-Sq analysis, V f = .60, AT = -500° F 



[ 0 ] 

[0 2 /±15] s 

[0 2 /±30] s 

[0 2 / + 45] s 

[0 2 /±60] s , 

[0 2 /±75] s , 

[0 2 /90 2 ] s 


Figure 42. Effects of laminate orientation on axial and radial Interfacial matrix stresses. 


P75/934, FE-Sq analysis, V f = .60, AT = -500° F 



Figure 43. Effects of laminate orientation on tangential and shear interfacial matrix stresses. 
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P75/934, FE-Sq analysis, V f = .60, AT = -500° F 
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[0 2 /±60] s 


x 9x 
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Figure 44. Effects of laminate orientation on out-of-plane shear interfacial matrix stresses. 


Table 1 1 shows a comparison of the maximum stresses predicted from the micromechanics FE 
analysis with those predicted from LPT for the [0 2 /902) s laminate. Both the maximum tensile and 
compressive values are shown, when present, for the FE predictions. LPT assumes that a state 
of plane stress exists in the laminate, with <r x . ay. and n<y as the only nonzero stress compo- 
nents. These stresses are computed for individual lamina, which are modeled as homogeneous 
orthotropic layers. The micromechanics FE analysis models the fiber and matrix as distinct indi- 
vidual phases, and the stresses predicted from this analysis refer to these individual phases. As 
shown in Table 11, the micromechanics FE analysis predicts a three-dimensional state of stress 
in the fiber and matrix, with stress levels much larger than the lamina components predicted from 
LPT. There differences have a significant influence on the prediction of thermally induced damage, 
as will be discussed in the next section. 


Results and Discussion 


71 



Table 11. Maximum stresses In [02/90 2 ] 8 P75/934 (^=0.60, /1T=-500°F ). 


LPT 


FE 

Lamina 

Fiber 

Matrix 

Stresses 

Stresses 

Stresses 

(Ksi) 

(Ksi) 

( Ksi) 

<7x — —9.95 

(t x = —30.3 

<r x = 19 9 

( 7y = 9.95 

<y r = 7.3, -8.9 

<r r = 121,-8.0 

Txy — 0.00 

<jq = 9.9, —4.5 
= 14, -7.6 

<rg = 23.4 
r rfl = 2.6, -9.3 


5.4 Thermally Induced Damage 

Thermally induced damage, in the form of matrix cracks, is well documented for polymer matrix 
composites exposed to temperatures far below their cure temperature^ 5-7,9,10 * . The thermal stress 
analysis presented in the preceding sections will be used to predict the initiation temperature 
and location of this type of damage. These predictions will be compared to experimental data 
generated for this study on unidirectional and [(V^ls P75/934 composites, and with literature 
data on T300/934 and T300/5208 composites. 

In addition to thermal stress field information, damage initiation predcitions also require failure 
strength data, and a failure criterion that relates the induced stresses to the failure strength. As 
stated earlier, the majority of experimental data suggests that thermally induced failures in polymer 
matrix composites take the form of matrix cracks. However, data on the failure strength of the neat 
matrix material is very limited. Elevated (+250 °F ) and RT tensile data on 934 epoxy are reported 
in the literature by Fox et. al. (59) , but no data for temperatures below RT could be found. The 
934 ultimate strengths reported by Fox et. al. were = 8 53 Ksi and = 10.30 Ksi at RT 
and +250°F , respectively. Fox et. al. also reported proportional limits, defined as the stress at 
which the stress-strain curve departs from linearity, of <r£] = 3.42 Ksi and <r™ = 2.44 Ksi at RT 
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and +250 °F , respectively. These data were obtained during static tension tests to failure (i.e. no 
unloading and reloading), and therefore it is not known whether this proportional limit corresponds 
to a material yield point. 

A maximum stress failure criterion, using the square array FE results, was used as the initial model 
to predict thermally induced failures in P75/934 laminates. This criterion predicts failure when 
any of the principal stresses first exceeds the ultimate strength of the material. For the case of 
failures in the 934 matrix, the maximum calculated tensile principal stress, ay , was compared to the 
ultimate tensile strength of the matrix, tx^. The principal stresses were calculated at temperatures 
of 250, 150, 75, -50, and -150°F which correspond to applied AT 's of -100, -200, -275, -400, and 
-500 °F , respectively. The results for a [0] and [02/902]s laminate at a V f = .50 are shown in 
Figure 45. The V* = .50 is consistent with the experimental P75/934 data (10) . The first failure in 
the [0] laminate was predicted at approximately 40°F . The maximum ay stress for this laminate 
always occured in the axial (x) direction on the fiber/matrix interface at the circumferential location 
shown on the figure. First failure in the [0 2 /90 2 ] s laminate was predicted at 95 °F . The direction 
of ay was in the transverse, or y direction, and again occured at the fiber/matrix interface at a 
circumferential location shown on the figure. Based on the results described in Section 5.3. 1.2 
and shown on Figure 30, failures in T300/934 composites would be predicted at approximately the 
same temperatures and locations. 

The results presented in Figure 45 do not agree with experimentally observed failures. First, 
no thermally induced failures have been observed in either P75/934, or T300/934 [0] laminates 
..^ected to repeated thermal cycling between ±250 °F < 4 ’ 10 ). Secondly, thermally induced failures 
in T300/934 and T300/5208 [0 2 /902] s laminates have only been observed in specimens exposed to 
temperatures below -200 °F (4 ’ 9) . T300/5208 is very similar in response and properties to T300/934. 
Thirdly, although P75/934 [O 2 / 9 O 2 ]® laminates do exhibit thermally induced failures at temperatures 
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P75/934, FE-Sq analysis, V f = .50 



Figure 45. Predicted thermally induced matrix failures using a maximum stress failure 
criterion. 


above RT* 10 \ the mode of failure appears to be due to a combination of a transverse component 
and a radial normal component at the interface. A photomicrograph of atypical thermally induced 
failure in a P75/934 [0 2 /90 2 ]s laminate is shown in Figure 46. 


Two possible explanations for the discrepencies between the experimental data and the results 
of Figure 45 are as follows. First, the FE stress analysis assumed linear elastic matrix behavior. 
However, as noted above, Fox et. al, reported a RT proportional limit of 3.42 Ksi in the stress-strain 
behavior of 934 epoxy. A von Mises yield criterion was used to predict the temperature at which 
this proportional limit would be reached. The von Mises criterion predicts yielding when 


r oct 


— \/2<7yp/3 


(5.4.1) 


where 


y/(g 1 - (T 2 )2 + (<t 2 - <T 3 ) 2 4- (<r 3 - <ri)^ 
r oct “ ~ ^ 


(5.4.2) 
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.002 h 

Aft«r fabrication 

Figure 46. Photomicrograph of typical thermally Induced failure In a P75/934 

laminate. 

and the components <r y , <r 2 , and a 3 are the principal stresses. The results predicted by substituting 
the RT value of <r£J for <r yp in equation (5.4.1) are shown in Figure 47. As shown in the figure, matrix 
yielding was predicted to begin at approximately 240 and 260 °F for the [0] and [ 02 / 902 ]* P75/934 
laminates, respectively. For the [0] laminate, the maximum value of occurred on the fiber/matrix 

interface at a circumferential location close to the y axis. For the [ 02 / 902 ]* laminate, the maximum 
value of root occured a short distance away from the interface at a circumferential location close to 
the 2 axis. The inclusion of this nonlinear matrix behavior in the analysis would have lowered the 
predicted stresses, and changed the predicted first failure loads. 

The second, and more conclusive explanation of the discrepancies between the experimentally 
obser thermal matrix failures, and the predicted failure loads and modes, is that the matrix 
ultimate strength, <rJJ, is not representative of the actual in-situ strength of the matrix. This claim is 
supported by the large differences observed between experimental values of < 7 ^, given above, and 
values of the transverse tensile strength (Yr) measured for unidirectional P75/934, at NASA-LaRC, 
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P75/934, FE-Sq analysis, V f = .50 



Figure 47. Predicted thermally induced matrix yielding using a von Mises yield criterion. 

and T300/934* 60 ). These values are given in Table 12. The difference in Y T for the two composite 
systems suggests a difference in the interfacial bond characteristics of the two composites. These 
data show that T300 fibers form a much stronger bond than P75 fibers with the 934 epoxy. The 
slightly larger value of Y T compared to at RT for the T300/934 composite is not understood, 
but may be within the scatter of the experimental data. One possible cause for this difference in 
Y t ’s is the difference in morphologies (Section 5.1) of these two fibers due to differences in their 
processing. Differences in Y T for composites with different modulus carbon fibers (AS4, 30 Msi 
fibers, and P75) have also been measured in thermoplastic PEEK matrix composites* 65 ). 


Table 12. Comparison of Neat Matrix and Transverse Lamina Ultimate Tensile Strengths. 


Temperature 

(°F) 

934 Epoxy 
(Ksi) 

P75/934 

Yj 

(Ksi) 

T300/934 

Y t 

(Ksi) 

250 

10.30 

3.20 

6.76 

75 

8.53 

3.50 

9.37 

-250 

NA 

1.00 

4.56 
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Predicted o r (FE-Sq analysis) 

---O Measured Yj 



Figure 48. Comparison of radial Interfacial matrix stress with measured transverse 
lamina strength. 


The maximum predicted matrix <r r interfacial stress as a function of temperature, for [0] and 
[0 2 /90 2 ] s P75/934 laminates, is compared to the experimental values of Y T in Figure 48. As 
discussed earlier, predicted <r r values for T300/934 and T300/5208 laminates would be very similar 
in magnitude to those presented in Figure 48 for the P75/934 laminates. Assuming Yj to be 
representative of the interfacial bond strength, failure (i.e. fiber/matrix disbonding) would be 
predicted when a, exceeded Y T . Based on this criterion and the results of Figure 48, failures 
were predicted in [0 2 /90 2 ] s P75/934 and T300/934 laminates at approximately 75°F and -1 50°F , 
respectively. No failures are predicted in [0] laminates of either system. These failure predictions 
are in excellent agreement with the experimental data cited earlier. 


The results presented in this section show that neat matrix strength is not a good measure of 
composite failure strength under thermal loading conditions. A failure criterion based on interfacial 
bond strength was much more successful in predicting thermally induced failures in the two 
composite systems studied. The lamina transverse strength, Yj, appears to be a good first 
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approximation of the interfacial bond strength. However, more research is needed to accurately 
quantify this strength. 
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6.0 Concluding Remarks 


The objective of this research study was to analyze the deformations and stresses, induced by 
thermal loading, in continuous fiber reinforced composites using the properties and behavior of the 
constitutents (i.e. fiber and matrix). The study focused on two primary areas. First, available explicit 
expressions for predicting the effective coefficients of thermal expansion (CTE's) for a composite 
were compared with each other, and with a finite element (FE) analysis that was developed as part 
of this study. Comparisons were made for a wide range of fiber/matrix systems to determine the 
influence of the constituent properties. All of the predictions were compared with experimental 
data to assess the validity and shortcomings of the individual analyses. The major conclusions 
from this portion of the study may be stated as follows. 

• All of the analyses investigated predicted nearly indentical values of the axial CTE, <*i , for a given 
material system, and all of the predictions were in good agreement with the experimental data. 
This implies that simple strength-of-materials (SOM) formulations are adequate for accurately 
predicting values of . 

• Results from the FE analysis, and those from the solution of a generalized plane strain boundary 
value problem, were in excellent agreement with each other, and with the experimental data for 
the transverse CTE, o 2 . The less rigorous plane stress and strength-of-materials formulations 
were in poor agreement with the experimental data for all of the material systems studied. 
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. As was expected, predicted values O «, were most sensitive to the axial properties of the fiber, 
and the predicted values o( „ were most sensitive ,0 matrix properties. the sensitivity ot o 2 to 
transverse fiber properties was less than expected. 

. Moderate amounts of cylindrical liber ortho, ropy (i.e. E, = 2E, or vice versa) had a negligible 
e«ec, on the predicted values o, e, and a. This degree of cy«ndricai Eber orthotropy is juried 
based on the excelled aggreem^t observed between predicted and expenmenta, values o, a, 
and o 2 when transversely isotropic liber behavior was assumed. 

The second portion ol this study focused on foe determination of thermally induced stress fields 
in foe individual constituents. Stresses predicted torn a FE anaiysis were compared to those 
predicted trom foe closed-form Composite Cylinder (CC) analysis for unidirectional laminates. The 
thermal load was representative of cooling from an elevated temperture cure (350- F ) to a typical 
cold (-tSO-F ) spacecraft environment. The constituent propenies and moosfotcturethat control 
foe behavior were identified and materia, modiflcatlons ,0 improve foe behavior were suggested. 

A Simpie global/local FE formulation was used to asses foe influence of loc* variations in fiber 
volume fraction, V-, and multiple ply laminate constraints on foe constituent stresses Thermal 
stress calculations were used to predict probable damage initiation locations, and the results were 
compared to experimentally observed damage in several epoxy mafox composites The maior 
conclusions from this portion of the study may be stated as follows 

. The type of analysis (i.e. CC or FE) significantly affected foe distributions and magnitudes of 
foe predicted thermally induced stresses. The most notable difference was foe absence of any 
tensile values of the radial matrix stress component, and foe absence ol any in-ptane shear 
stresses. These differences are due to foe fact that the CC analysis does not account for the 

influence of adjacent fibers. 
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Thermally Induced marix stresses were no, a strong lunclion o, fiber properties. Conseco 
moderate amounts of fiber orthotopy had negligible effects on the malrix stresses. The two most 
promising rnemods ,o «««* 'owe, the magndudes of themra* Induced madx stresses 
are reducing the product of die matrix modulus and CTC (E M n m ). and reducing the applied tem- 
perature Change. For a .bred opera, ng temperaae. this Implies lowering the cure temperature. 

. Manx aresses Increased in absolute value with increasing V. Globaldoca, ana*ses showed 
tha aresses in local regions of high V. constrained by a lower global V, can be approximated 

by an unconstrained analysis of the high V f region. 

. Multidirectional [0*/ * «b 'aminaes had larger predicted marix aresses than unidirectional ([0]) 
laminaes. The stresses increased with increasing lammaion angle 0, up to 
magnitudes o, matrix aresses remaned approximaay consUm, ,o, * > 60- up to „ = 90- . The 
most significant deferences from the [0] laminae response were the presence of a lage tensile 
radial marix stress component, a, , a, the intedace. and the exiaence of non-zero out-of-ptane 

matrix shear stresses, r ex and r xr . 

. Thermaiy induced matrix lailure predictions, using a maximum stress Mure criterion and die 
ultimae tensile arength o, the malrix, „ 5 . were in very poo, agreement with experimenta daa. 
The wo reasons poaulated for these discrepancies were non-linea marix behavior, and <r& a 
nea matrix propeny. no. being represent of the in-situ marix strength. The second reason 
was supported b, experimenta, obsecrations of marix damage, and measured vaues o, trans- 
verse lamina strength. Y T . A lailure erdehon based on fh. radia, intedacia stress component 
and Y t , was in excellent agreement with experimental data. 

summary, a micromechanics anaysis was devaoped. and used for predicting the therma* in- 
duced defomtaions and aresses in a composite exposed to a typica spacecraft environment 
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Predicted deformations of the composite were in good agreement with experimental data. Cal- 
culated matrix stresses were used to predict thermally induced matrix failures, and were in good 
agreement with the experimental data when representative values of the in-situ matrix strength 
were used in the failure criterion. Future research should focus on more accurate modeling of 
matrix behavior (i.e. nonlinear behavior), and a more thorough understanding of interfacial bond 
characteristics. 
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Appendix A. Rosen and Hashin Analysis 


Due to the transverse isotropy of the constituents and the unidirectional composite, the only non- 
zero components in the material property tensors of equations (2.1.8) and (2.1.9) are given by* 23 ) 


ft 11 = 0 ( 1 ), «22 = «33 = 0 ( 2 ) 

= 1/E ( i), S2222 = S3333 = 1/E(2, 

S 1122 = S1133 = -i/(,2)/ E (1)l S2233 = -"( 23 )/ E ( 2 ) (A. 1 ) 

S1212 = S1313 = 1 / 4 G(, 2 ), 3 2233 = 1 / 4 G ( 23 ) 

These relations are true for the effective property (superscript ' ), the volume average property 
(superscript ), and the constituent property (superscripts f.m) tensors. The subscripts on the 
engineering constants refer to the material property coordinates of Figure 4, and are enclosed in 
parentheses to distinguish them from tensor indices. 

Using contracted tensor notation* 66 ) (subscripts a,b = 1,2, 3, 4,5, 6), equation (A.1) may be rewritten 
as 

“I = 0(1), r >2 = o 3 = ft<2) 

Appendix A. Rosen and Hashin Analysis 8Q 



S u = 1/E (1 ), S 22 = s 33 - V E (2) 


Si2 =S 13 = -W E (1)' S 23 = -^(23)/E(2) 

S 66 = S 5 5=1/4G(12), S 4 4 = 1/4G (2 3) 


Now, equations (2.1.11a) and (2.1.11b) may 


be rewritten using contracted notation as 


«i 


= a, + (a^ ~ a^ m) )P a b(Sb1 ~ § bl) 


(A. 3) 


«2 


= « 2 + K' - 


<? - ai m) )Pab(Sb2 - Sb2) 


(A -4) 


Equations (A.3) and (A.4, may be expanded and singed by ^nq on* the nonzero terms and 
noting that S,„ is symmetric and theretore P* is symmeaic. Tbe Wowing relations can also be 

shown for P a b 


P 22 = P33, Pi2 = p i3- Pss-Pee 
Equations (A.3) and (A.4) may now be written as 


(A. 5) 


t*1 


+ (S„ - SnlHnf 1 - efVii + ("2 

+ (S, 2 - $12)1(0? - of )2 p 12 + <"2 


- 4 m) )2P 12 ] 

- af" , )2(P 22 + p 23)) 


(A 6) 
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+ <§* - Sa)t(«? - «rv« + <«? - "2 H )C*» + P *» 

+ A.-^f-r^«^- -»?•><*. +«*» 


(A 7) 


. /o i q\ which may be rewritten in contracted form as 
where P,„ is determined from equal, on (2.1 .9) *»<» ™» 



(A 8) 


•nte solution ot equation (A.8) gives the elements ot P* 
and may be written as 


appearing in 


equations (A. 6) and (A. 7), 


P 11 = (A 22 A 2 2-A 23 A23)/DetA 

P 12 = (A 12 A 23 - A22Ai 2 )/Det A 
p 22 = (An A22 " Ai 2 Ai 2 )/ Det A 
p 23 = (A 12 A 12 - An A23)/Det A 


• ^ r„nH a c _ s'T'), and Det A is given by 
where A ab is defined as (b ab o ab 


(A. 9) 


Det A = lA^A* - A 23 A 23 > + 2A, 2 (A, 2 A 2 3 - A*A, 2 H 


(A. 10) 
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Appendix A. Rosen 


and Hashin Analysis 



Appendix B. Finite Element Formulation 


As stated in Section 3.4, 8-node quadratic isoparametric elements were used for the FE analysis. 
The global and local coordinate systems for this element are shown in Figure B1. The relationship 
between the global and local coordinates can be written as 

8 

y = N iV< (B 1 ) 

i-1 

8 

z = X] N i 2 i (B.2) 

i-t 

where the interpolation polynomials N| are functions of the local nodal coordinates ((, v ), and yj 
and Zj are the global coordinates of the i’th node. 

The interpolation polynomials for the 8-node elements are quadratic, and defined as 


N i - 4O + &)(1 + r/ 0 )(so + 7/0 - 1), 

i = 1,3,5, 7 

(B.3) 

N i = ~ ^ 2 )( 1 + Vo), 6 = 0 , 

i = 2,6 

(B.4) 

Ni = ^( 1 +eo)( 1 - 7 / 2 ), r H = 0, 

i = 4, 8 

(B.5) 
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y 

Global coordinates (x, y, z) 
Local coordinates ( ^ ) 


x 

Figure B1. Element local and global coordinate systems. 


where rj 0 = tjrji, and £ and 7 range from -1 to +1 at the corners. 


Isoparametric elements use the same interpolation functions for the geometry and the nodal vari- 
ables (i.e. displacements). Therefore, the unknown displacement functions of equation (3.4.1) 
may be written as 

8 

U(y,z) = £NiUi 


a 

V(y,z) = £NiVj (B.6) 

i=1 

6 

W(y,z) = ^NjWj 

i= 1 

where u h v jf and Wj are the unknown displacements at the nodes. 


A displacement based FE formulation requires the minimization of the total potential energy, PE, for 
each element with respect to the unknown variables. This procedure leads to a set of linear simul- 
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taneous equations relating the unknown nodal displacements to the applied loads. The equations 
tor each element are assembled into a global system, the boundary conditions are applied, and 
me system is solved tor the nodal displacements. A detailed description of this procedure may 
be found in the literature^. An abbreviated description ol this general procedure, and special 
modifications for the specific problem under consideration, follows. 

The total PE of an element is the sum of the strain energy, U s , and the work of the externa 
W L . The strain energy for an element is written as 

U s = 1 / W{f M }dvol (B7) 

JVOl 

where M and {<„) are the element stresses and mechanical strains, respectively, and have the 
form 

T \ (B.8) 

{<t} T = {<X X , <Ty , <T Z , Tyz- r ZX, T *Y » 

{<m) T = {fx.fy> f z,7yz,7zx.'bcy} 


The work of the external loads is written as 

W L = — {q} |f } T — f x^x 


(BIO) 


where {q} are the nodal displacements, {f} are the applied nodal loads, and F° is the uniform 
average axial force (i.e. a scalar quantity). 


The total strain for a linear elastic system is written as the sum of the mechanical 


strain and the free 


thermal strain, or 


{(} = {<m) + {<t} 


(B 11) 


Rearranging to solve for the mechanical strain leads to 

(fM) = M - ( ( t) 


(B 12) 
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Minimization of the total PE with respect to the nodal displacements requires that {<r} and {f m ) 
in equation (B.7) be expressed in terms of the displacements, and is accomplished by using the 
strain-displacement relationships of linear strain theory. These relationships are written in matrix 
form as 

w = nm (®.13) 

where {u} is the vector of total displacements, and the matrix [L] is defined as 


rd/dx 

0 

0 

0 

d/dy 

0 

0 

0 

d/dz 

0 

d/dz 

0/d y 

d/dz 

0 

d/dx 

lO/dy 

d/dx 

0 


Note that engineering shear strain has been used in the above relationships. 


Combining equations (3.4. 1 ) and (B.6) and substituting into equation (B.1 3) leads to an expression 
for strain in terms of the nodal displacements written as 

{0 = [B] {q} + {£} (B 15) 


where 

{q} T - {u 1 ,v 1 ,w 1 ,...u 8 ,v 8 ,w 8 } 


and 


K) T = {•?. 0.0, 0,0,0} 


(B 16) 


(B 17) 


The matrix [B] (6 x 24) is defined as the matrix product of [L] and a partitioned matrix involving the 
interpolation polynomials (equations (B.3- B.4)) written as 

(B) = «L]) ([l](N 1 )|[l](N 2 )| . . . [I] (Nq)) (B 18) 
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where I is a 3x3 identity matrix. Substitution of equation (B.15) into equation (B.12) leads to an 
expression for the mechanical strains in terms of the displacements, written as 

+ < B19) 

where the free thermal strain vector is defined as 

{f T } = {a}zlT ( B - 2 °) 


and 


{a}' = {«x> a y> 0 'z> Q 'yz> a zx> Q xy} 


(B.21) 


and 4T is a uniform change in tempertaure. 


The stresses are expressed in terms of the displacements through the material constitutive equa- 
tions, written in general form as 

M = [C]{r M } (B.22) 


The specific form of [C] depends upon the relative orientation of the material principal axes with 
respect to the global (x,y,z) axes. The elements of [C] for an orthotropic material (i.e. a material with 
3 principal planes of symmetry) with two different orientations (coincident , and a rotation about the 
x-axis) are given in Appendix C. For rotations about the x-axis, the angle of rotation, 6 , varies with 
position in the yz plane, and therefore varies within a given element. This variation is included in 
the analysis by computing 0 from the global coordinates at each of the Gauss points used for the 
numerical integration in equations (B.26-B.27) and equations (B.30-B.32). The matrix {a}, defined 
in equation (B.21), also depends upon the relative orientation of the principal material axes with 
respect to the global axes. The elements of {»} for the two orientations described above are also 
given in Appendix C. 
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Equations (B. 1 9) and (B.22) are substituted into equation (B.7) to write an expression for the strain 
energy in terms of the displacements as 


U s = l [ {8} T [C]{B}dvol 
d Jv ol 


(B.23) 


where 


{B} = ([B]{q} + K}-{ fT }) 


(B 24) 


Two cases must be considered for the minimization of the total PE (equations (B.10) and (B.23)) 
with respect to the unknowns. First, for the case of a known t £, the PE is minimized with respect 
to only the nodal displacements, {q}. This minimization results in an expression written as 

[K]{q} = {F} (B.25) 


where 

[K] = / [B] T [C][B]dvol 

Jv ol 


(B.26) 


and 


{F} = - / [B] T [C]({<®} - {f T })dvol + {f} 
Jv ol 


(B 27) 


The case of an unknown (° requires that the PE be minimized with respect to both the nodal 
displacements, {q} , and This formulation is specific to the case of generalized plane strain with 

an applied uniform axial force, F°, and does not appear in the reference* 41 ) cited previously. The 
minimization of the PE for this case results in an expression of the general form 


[K*](q* } = {F-} 


which may be expanded and written as 




(B.28) 


(B.29) 
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when, |K], {q>, {F>. and F x ° have the same definitions as previously given, end 

{K x } = 1'st column of j [B] T |C]rhol 

K°= [ [CnJdvol 

Jvol 


and 


F x = 1'st element 


of / [C]{< 

e/VOl 


T }dvol 


(B.30) 

(B.31) 

(B.32) 


Two steps are requiredto compute theelementsof (Kl, (K*l, {F),and{F*}. First, the derivatives ol 

me interpolation polynomials, N„ with respect to the global coordinates, (y.z). must be evaluated 
lor terms involving [B|. Howler. the shape .unions are in terms olthe lodgement coordint.es, 
rivative evaluation is acomplished by use ol the Jacobian matrix in an expression 


(£>»?)• The derivative 
written as 


I OKJdz J 

where the Jacobian matrix is defined as 

[J] = 




(B.33) 



. 3N 8 /^1 

-yi z r 


. dN 0 /d»?j 

-ya Za- 


(B.34) 


The second step is to numerically evaluate the integrals in equations (8,26- B.32). This is accom- 
plished by first translorming the integral into the (£ , ,) coordinate system noting that 

dvol = |J|(d£)(dr/) (B.35) 

where |J| is the determinate ol the Jacobian matrix. The limits of integration become -1 to +1 on 
bom e and Gauss quadrature is used lor the integration and all ol the integrals can be written in 


the form 



n n 


Gdvol^^^HiHjG^.r^lJI 

i=1 j=1 


(B.36) 
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, a(f , Tl .\ j S the funciion to be integrated. A 2x2 (n 2) 
^ u anH H- are the weight functions end G(£i> 

zzz — - — 


Gauss points are given in the 


FE reference text cited earlier* 


Equations (B.25) and (B.29) represent a system ot linear 


simultaneous equations for one element. 


element equations are assembled into a 


global system ol equations by requiring continuity 


These The prescribed boundary conditions 

cme displacements a, the nodes ol adding elements. pr 
are imposed and the resulting system of simultaneous linear equations is solved for the 
nopal dispalcements. are. it -essrey «J. Tbe element s.r*ns and stress are comp 
the nodal displacements by using equations (B.13) and (B.22). 

Resultant internal nodal forces are computed from 

, _ , (B.37) 

{ P) = / [Bl T {<r}dvol 

Jv ol 

Where P is the vector of resultant nodal forces for an element. The integration in equal, on (B.37) rs 

evaluated numerically by rewriting Sin tbe form of equation (B.35). The summation of theseno 

forces in the x, y. and z directions, respectively, must equal zero for static equilibrium, and can be 
used as a check on the validity of the solution. 

tion presented above (U. ntiremlzation of the total potential 


It should be noted that the formulation 


that [K] and [K*l will be symmetric 


and positive definite, after the imposition of 


exists and a solution can be found. 

The three types ot boundary conditions allowed in the displacement formulation presented ab0 *® 
are prescribed nodal displacements, applied nodal forces, and constrained nodal 
The firs, hvo involve steward „ocedures fo, approptiar* modifying the global - 
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(IK] or IK*]) and load vector ({F} or {F*}), and are well documented in texts on finite element 
analysis! 41 ). The application of constrained nodal displacements is not as well documented and 
requires further explanation. Constrained nodal displacements refer to a series of nodes whose 
displacements are proportional to one another but whose absolute values are unknown (e.g. v, = 
v 3 = 4v 6 ). This effectively reduces the total number of unknowns in the system. The global stiffness 
matrix and load vector are modified by combining the contributions of the affected nodes into a 
single row and/or column. These modifications do not destroy the symmetry of the stiffness matrix, 
but can increase the bandwidth (to be discussed subsequently). The specifics of the procedure 
for modifying the stiffness matrix and load vector are given by Adams and Crane (26 >. 

The solution of equations (B.25) or (B.29) involves solving a system of simultaneous linear alge- 
braic equations. As stated earlier, [K] and [K* ] will be symmetric and positive definite, and therefore 
an inverse will exist, and a solution is possible. The matrix [K] will also be banded due to the as- 
sembly procedure in which a nodal displacement affects only those elements adjoining that node. 
However, [K*] is not banded due to the fact that every element contributes to f °. This type of ma- 
trix is sometimes refered to as an ’ arrowhead" matirx. Both matrices are shown schematically in 
Figure B2. Numerous solution algorithms exist for banded symmetric systems. These algorithms 
require storage of only the terms that lie on or above the diagonal and within the band, thus greatly 
reducing the amount of computer memory and time needed for the solution. A special solution 
algorithm was developed for this study that allows a banded storage scheme to be used for [K*] 
also. This was accomplished by storing the terms of [K* ] which contribute to^asa separate vec- 
tor. A standard [L] [U] decomposition was then used to solve the banded system for the unknowns, 
with terms from this separate vector used when needed. 
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Figure B2. Finite element global stiffness matrix architecture. 
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Appendix C. Finite Element Constitutive Equations 


Ttre material sMnoss matrix, [d is a 4 '.h order symmetric tensor. Fot, owing the contracted tensor 
notation used in Appendix A. Id can be wriden as a 6x6 symmetric matrix as 


Cn C 12 C 13 0 0 0 1 

C 12 C22 C 23 0 ° 

C,3 C 23 C 33 0 0 0 

o 0 o C44 0 0 

0 0 0 0 C55 0 1 

0 0 0 0 0 C66J 


(C. 1 ) 


where the elements of [C] are given by 


1-^32 = !ji2 + ^32^13 c 13 = ^^ 1 7r 1 

°ii - - e 2 E 3 X~ ° 12 E,E 3 4 

1 - ^13^1 c 23 =^^ C 33 = 1 i^ff i 

c 22 - EtE 3 zi 23 ^1^2^ 1 2 

C44 = G 23 C55 = G3I c 66 = °12 


(C. 2 ) 


and 


1 _ - V2 3^32 ~ ^ 13^31 

A = Ei E 2 E 3 


(C. 3 ) 


The terms Ej. Gy, and ^ in the above equations 
Poisson’s ratios, respectively, in the material principal coordinates. 


refer to the Young’s moduli, shear moduli, and 
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material principal (1,2,3) axes are coincident 


with the global (x.y.z) axes, then the [C] 


most Ol me micromechanics problems analyzed ,n m,s study. 

However materials that possess cylindrical orthotropy (circumferential or radial) in the (2,3) 
require a transformation of the IC] matrix to (Cl- The appropriate trarwtormatiornnvolves 
about me global z *ls, and Is s^wn - Hgure C, . Tbe elements o, ,0, are 

mined using a standard 4‘th order tensor transformation law t66) . and can be written in contra 


notation as 


Cii = 


C„ = Ci i 

C 12 = m 2 C 12 + ° 2 ^i3 
Ci3 = n 2 C 12 + m 2 C 13 
C 14 = mn(C 12 -C 13 ) 

Cjj = m'Cjj + mVpCja + 4C44) + P <c 33 
c 23 - mVfC* - 4C« + 0 M ) + (rn 4 + n 4 )^ 

C 24 = m 3 n(C 22 - C 23 - 2C„) + * 2 C»-C X ) 

C M > m 4 C 3 3 + + 4C "> + n c zz 

C54 - m 3 n(Ca + 2C„ - C M ) a- ™ ! fe ' °*> “ 2C “ ) 
C„ = m 4 C 44 + m*n?(C* - ^ - 2C„ + C33) + n 4 C 44 

C55 = m 2 C55 + n 2 C 6 6 

C56 = mn(C66 - C 55 ) 

Cqq = m 2 C 6 $ + n 2 C 55 


where all other C, - 0, and m = cos» and n = sin«. 
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Figure Cl . Cylindrical material transformation geometry. 

The elements of {«} used in Appendix B, for orientations where the material 
coincident with the global axes, are defined as 

«x = «i , a y = o 2 , a z = « 3) 

«yz = 0, »zx = o, ftxy = o 

where ai , « 2 , and « 3 are the CTE’s in the material principal directions. 

The elements of {<*} for a rotation about the global x axis are given by 

<*X = <*i 

Gy = m^(*2 + n 2 03 
«z = n 2 r* 2 + m 2 o 3 
Gyz — 2mn (a 2 — 03) 

«zx = 0 
«xy = 0 

Appendix C. Finite Element Constitutive Equations 


principal axes are 


(C5) 


(C.6) 
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I 


■nw above transformation for (n) is based on the tee of engineering shear strain in the FE formu 
lation, as opposed to tensor shear strain. 
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Appendix D. Composite Cylinder Solution 


The orthotropic constitutive relations in terms of the cylindrical coordinate system of Figure 7 may 
be written as 


( cr x 


fCn 

Ci 2 

Cl 3 

0 

0 

0 - 


' e x — AT ' 



Ci 2 

C22 

C23 

0 

0 

0 


€> OC2AT 

J <T$ 

y __ 

c 13 

C 23 

C33 

0 

0 

0 

J 

£9 — a$AT 

j r re 


0 

0 

0 

C44 

0 

0 

\ 

Ire 

I T »x 


0 

0 

0 

0 

C55 

0 


7BX 

l r X r > 


. 0 

0 

0 

0 

0 

C66- 


[ 7xr « 


(D 1) 


where the [C] matrix has the same definition as given in Appendix C, except that the material 

principal coordinates 1 ,2, and 3 are coincident with the x, r, and 6 directions, respectively, of Figure 
7. 


The solution to equation (3.3.2.S), presented by Avery and Herakovich* 35 ). has two forms depend- 
ing upon whether the material is cylindrical^ orthotropic (i.e. properties in r and 6 directions are 
not equal), or transversely isotropic (i.e. properties in r and 6 directions are equal). If the fiber is 
cylindrical^ orthotropic the displacements in the fiber are given by 


where 


w f (r) = A' r A ' + A^r A z + LjeJr + L^T 



(D.2) 


(D.3) 
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and 


i f _ (Qr — ^'to) Q 'i 

4 - C^-C M (0.4) 


II the liber is transversely isotropic the displacements in the liber are given by a much simpler 
equation written as 

w , (r) = A{r V . +A»r A 2 ( D . 5) 

The term A, 2 , in both of the above displacement fields, is given by 



The displacements in the matrix, assuming isotropic material behavior, take the same form as 
equation (D.5) and are written as 


w m (r) = A'V7 , +A£V? (D7) 

where A 1 2 is given by equation (D.6) using the stiffness coefficients of the matrix phase. For 
isotropic or transversely isotropic phases the values of A are given by A, 2 = ±1. 

The stresses in each phase are obtained by substituting the the expressions for the radial dis- 
placemts, equations (D.2,D.5,D.7), into the first equilibrium equation of Section 3.3.2 (equation 

(3.3.2.3)), and the results into equation (D.1). The resulting expressions for the stresses in each 
phase are given by 


'i = A* (Cl + c£A f ,)r A «-i + A»(C{, - C' r A')r^-’ + M \ ( - + N [_iT (D.8) 


with 


M| = C? x + L'(C[ tf+ c! r ) 


(D 9) 
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and 


n| = 4(cI+ c I)- c V v ! 


(DIO) 


for a cylindrical^ orthotropic fiber, and 


a\ = <Cj 9 + C- r ) + A<(Cj - C? r ) ~ + C' K ,° - C\a\AT 


-f f . 


(Dll) 


for a transversely isotropic fiber, and 


ff [" = A™(C£ + C") + A^(C£ - C-)i + C^° - 


(D 12) 


for an isotropic matrix. 

.. c A f Am A m a nd <? are determined from the following 
The five constants in the above equations, A, , Aj, i > 2 • 

boundary and continuity conditions. 

r _ 0. This in conjunction with the fact that A f 2 < 0 


1. The radial displacement w must be zero at 1 
requires that Ai, be equal to zero. 


! across the fiber/matrix interface, w f (a) - w m (a). 


2. The radial displacement must be continuous ; 

3. The radial component of stress, <r r , must be continuous across the fiber matrix interf 
alia) = a r m (a). 


4. There are no tractions applied on the outer 


therfore 


boundary of the matrix for pure thermal loading, 
= 0. The condition that r™ (b) = 0 *as already used .0 determine .he constant 
of the second equilibrium equation (Section 3.3.2). 


of integration arising from the solution 

5. The net axial force on the composite must be zero for pure thermal loading, or 
2 tt J 0 b <r x rdr = 0. 
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Conditions 2-4 result in 4 equations that can be soived sinuitaneeus* lor the remaining unkowns. 
A*, A™, AIJ 1 , and c?. 
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